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Abstract 

Air-to-air  missile  endgame  programs  typically  require 
tens  of  seconds  of  computer  time  to  calculate  a  single  esti¬ 
mate  of  probability  of  kill  (Pj,).  Thousands  of  these  esti¬ 
mates  may  be  required  in  any  given  endgame  analysis.  The  main 
objective  of  this  research  was  to  develop  a  method  for 
characterizing  an  endgame  model  using  only  a  fraction  of  the 
amount  of  data  generated  by  such  analyses.  In  pursuit  of  this 
objective,  a  second  objective  of  actually  accomplishing  a 
characterization  of  an  existing  endgame  program  was  established. 

The  endgame  program  SHAZAM  was  used  to  generate  data 
upon  which  to  base  a  characterization.  A  purely  empirical 
approach  was  first  attempted  wherein  all  main  and  interactive 
effects  of  the  critical  factors  were  to  be  assessed  strictly 
by  experimentation  with  and  observation  of  the  behavior  of 
SHAZAM.  The  high  order  of  interactions  between  the  factors 
precluded  completion  of  this  approach. 

A  semi-empirical  approach  was  then  pursued  whereby  the 
interactions  between  the  endgame  geometry  and  warhead  charac¬ 
teristics  were  accounted  for  geometrically.  A  two-parameter 
exponential  function  was  used  to  describe  the  behavior  of  Pj, 
as  a  function  of  the  warhead  burst  point.  Values  of  the  two 
parameters  were  calculated  based  upon  the  settings  of  the 
remaining  factors.  A  "vulnerable  length"  was  used  to  account 

xi  i 


for  the  target's  vulnerability.  After  several  iterations  of 
refinement,  the  characterization  was  compared  to  5HAZAM  for 
ten  randomly  selected  trajectories.  The  characterization's 
Pj,  estimates  were  within  0.1  of  SHAZAM's  estimates  for  97  per¬ 
cent  of  these  random  data.  A  simple  computer  program  imple¬ 
menting  the  mathematics  of  the  characterization  required  a 
maximum  of  0.2  seconds  of  computer  time  per  estimate. 
Approximately  two  hours  of  computer  time  were  consumed  by 
SHAZAM  in  generating  the  necessary  data.  Thus,  with  a  rela¬ 
tively  modest  investment  of  computer  time,  two  orders  of 
magnitude  improvement  in  time  required  per  P„  estimate  was 
achieved. 

The  suggested  method  for  characterizing  air-to-air 
missile  endgame  programs  is  based  on  the  semi-empirical 
approach  and  its  application  to  SHAZAM.  The  method  is 
heuristic,  rather  than  algorithmic,  due  to  the  subjective 
nature  of  the  characterization  process. 


A  METHOD  OF  CHARACTERIZING  AIR-TO-AIR 
MISSILE  ENDGAME  MODELS 

I .  Introduction 

Background 

The  last  100  milliseconds  or  so  of  an  air-to-air 
missile's  encounter  with  a  target  are  critical  to  its  even¬ 
tual  success  or  failure.  The  guidance  system  has  already 
done  what  it  is  capable  of  doing  and  can  no  longer  change 
the  outcome  of  the  encounter.  If  it  has  done  its  job  well, 
the  missile  will  hit  the  target  and  will  be  considered 
successful.  More  likely,  at  least  for  most  missiles,  is  that 
the  missile  itself  will  not  make  physical  contact  with  the 
target.  In  this  case,  it  is  up  to  the  fuze  and  warhead 
system  to  destroy,  or  "kill,"  the  target. 

The  warhead  is  a  combination  of  high  explosive  and 
metal  and  can  kill  the  target  with  one  of  two  mechanisms. 

The  first  is  with  the  air  blast  accompanying  warhead  detona¬ 
tion.  If  detonation  occurs  sufficiently  close  to  the  target, 
tli e  resulting  overpressure  will  cause  structural  failure  and 
subsequent  destruction.  Otherwise,  destruction  can  only 
result  from  fragmentation,  the  second  kill  mechanism.  Upon 
detonation,  the  metal  in  the  warhead  is  fragmented  and  ejected 
at  very  high  velocities.  If  one  of  these  fragments  strikes  a 
vulnerable  component  of  the  target,  there  is  a  finite 
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probability  that  the  entire  target  will  be  destroyed. 

The  purpose  of  the  fuze  is  to  detonate  the  warhead  at 
a  time  which  allows  at  least  one  of  the  warhead  kill  mechan¬ 
isms  to  destroy  the  target. 

Analysis  Methodology 

Over  the  years,  much  analysis  of  the  terminal  encounter, 
or  "endgame,"^  situation  has  been  done.  Numerous  computer 
programs  have  been  written  which  predict  missile  kill  proba¬ 
bilities-  -ENDGAME ,  AMEGS ,  SCAN,  SHAZAM ,  and  SESTEM  II,  just 
to  name  a  few  (Ref  1:  Sec  5-40  to  5-55).  Not  all  of  the  pro¬ 
grams  deal  with  the  endgame  problem  at  the  same  level  of  detail, 
and  some  of  them  perform  other  functions  as  well. 

There  are  numerous  purposes  for  which  endgame  analyses 
are  undertaken.  One  might  be  trying  to  select  the  best  war¬ 
head  out  of  a  number  of  candidates,  or  he  might  be  trying  to 
find  ways  to  improve  the  effectiveness  of  an  existing  one. 

Other  objectives  might  be  estimation  of  an  existing  fuze/ 
warhead  system's  effectiveness,  generation  of  effectiveness 
tables,  or  selection  of  a  fuze  time  delay. 

Regardless  of  the  objective,  the  underlying  question 
being  asked  in  any  such  analysis  is  this:  given  a  specific 
terminal  encounter  situation,  what  is  the  probability  that 
the  target  is  killed  by  the  missile?  This  question  may  be 

■'"Henceforth ,  the  term  "endgame"  is  used  to  mean  that 
part  of  the  missile's  trajectory  during  which  the  closest 
approach  to  the  target  occurs  and  all  velocity  vectors  can 
be  considered  constant. 
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asked  thousands  of  times  during  a  given  analysis.  Before 
the  question  can  be  answered,  a  more  fundamental  question 
must  be  answered:  what  are  the  vulnerabilities  of  the 
target?  For  any  given  target,  the  answer  to  this  question 
is  very  difficult  to  obtain.  Usually,  teams  of  vulnera¬ 
bility  analysts  combine  test  data  with  their  knowledge  of 
physics,  mechanics,  material  properties,  and  hydrodynamics 
to  create  what  is  known  as  a  "vulnerability  model."  These 
"models"  come  in  a  variety  of  forms,  and  generally  an  endgame 
program  is  structured  to  accept  only  a  limited  number  of  these 
forms . 

The  purpose  of  an  endgame  program,  then,  is  to  account 
for  the  interaction  between  the  warhead,  the  dynamics  of  the 
terminal  encounter,  and  the  vulnerability  model  to  yield  an 
estimate  of  the  probability  of  kill  (P^)  . 

The  most  common  approach  to  modeling  the  endgame  situ¬ 
ation  is  to  describe  the  physical  process  of  fragment  ejection 
mathematically.  Two  binary  criteria  are  examined  first: 
whether  or  not  the  encounter  was  a  direct  hit,  or  whether  or 
not  detonation  occurred  within  the  blast  envelope.  If  neither 
criterion  is  met ,  then  the  dynamics  of  the  encounter  are  com¬ 
bined  with  the  motion  of  the  fragments  to  determine  where  and 
how  often  the  fragments  strike  the  target.  Frequently, 
there  are  stochastic  elements  in  the  mathematical  description 
of  fragment  ejection  and  motion.  Also,  fragmentation  may  be 
modeled  as  a  "beam"  with  a  certain  density  of  fragments.  In 
this  case,  the  area  of  the  intersection  of  the  beam  with  the 


target  is  calculated.  In  either  case,  the  frequency,  loca¬ 
tion,  and  characteristics  of  the  fragments  striking  the  target 
are  compared  to  or  combined  with  the  vulnerability  model  to 
yield  an  estimate  of  P^.. 

Frequently,  endgame  programs  are  operated  in  a  "Monte 
Carlo"  fashion  in  which  variables  such  as  miss  distance  are 
treated  as  random  variables  and  sampled  from  distributions. 

Each  iteration  results  in  a  kill  probability.  The  average 
of  these  individual  probabilities  is  the  overall  probability 
of  kill  for  the  encounter.  The  sampling  variance  of  the  mean 
is  reduced  by  increasing  the  number  of  iterations. 

One  way  of  viewing  an  endgame  program  is  that  it  is 
just  a  "black  box"  (see  Figure  1).  Certain  inputs  are 
required,  the  program  processes  the  information,  and  a 
comes  out.  Most  endgame  "black  boxes"  have  two  character¬ 
istics  in  common:  they  are  "large"  in  terms  of  computer 
memory  required,  and  they  are  slow  (and  consequently  expensive 
to  operate)  .  Ten  seconds  of  central  processor  time  per 
estimation  is  considered  "fast"  in  the  endgame  business. 

Problem  Statement 

Given  that  an  endgame  program  is  run  thousands  of  times 
per  analysis,  is  it  possible  to  characterize  the  behavior  of 
the  output,  Pj. ,  with  respect  to  changes  in  the  inputs  so  that 
this  characterization  can  be  used  for  subsequent  P^  estima¬ 
tions?  If  so,  how  does  one  develop  such  a  characterization? 

It  is  the  purpose  of  this  thesis  to  investigate  these  questions. 
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Scope 


To  answer  the  above  questions  for  the  general  case  is 
a  large  and,  perhaps,  impossible  task.  This  thesis  investi¬ 
gates  the  question  for  the  case  of  a  single  endgame  program 
with  a  single  vulnerability  model  and  a  single  warhead 
description.  The  role  of  the  fuze  is  not  considered  in  the 
investigation.  A  characterization  accomplished  under  these 
circumstances  would  indicate  that  it  is  possible  the  same 
method  can  be  applied  to  different  programs,  vulnerability 
models  and  warhead  descriptions, 

Ob j  ectives 

Given  the  nature  and  scope  of  the  problem,  two  specific 
objectives  must  be  accomplished.  One  is  to  obtain  a  "charac¬ 
terization"  of  an  endgame  program's  behavior.  The  other  is 
to  define  a  method  for  arriving  at  this  characterization. 

By  characterization  is  meant  an  entity  which  mimics 
the  behavior  of  the  endgame  program.  There  are  many  forms 
that  such  an  entity  can  have:  equations,  statistical  rela¬ 
tionships,  algorithms,  and  tables,  just  to  mention  a  few. 

For  purposes  of  this  thesis,  the  tabular  ±orm--i.e.,  tables 
of  outputs  versus  inputs- -is  unacceptable  because  it  would 
not  be  an  improvement  over  existing  methodology.  Most  any 
other  form  should  be  acceptable,  as  long  as  it  is  not  more 
cumbersome,  expensive  or  time-consuming  than  the  "black  box" 
it  replaces. 

By  mimic  is  meant  that  given  the  same  inputs,  the 
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characterization  should  yield  the  same  output  as  the  endgame 
program  would  have.  Exact  duplication  of  results,  however, 
is  a  rather  heroic  objective.  A  more  realistic  goal  is  to 
yield  an  output,  which  is  within  an  acceptable  tolerance  of 
the  real  output.  As  mentioned  previously,  the  output  of  many 
endgame  programs  has  a  statistical  uncertainty  associated  with 
it  which  is  dependent  upon  the  number  of  iterations  run.  This 
is  not  an  uncertainty  about  the  real  value  cf  ,  but  is  an 
uncertainty  aboui  the  F„  the  program  would  predict  with  a 
very  large  number  of  iterations.  This  uncertainty  is  usually 
expressed  in  the  form  of  a  confidence  interval  about  the  Pj, 
estimate.  A  commonly  accepted  level  of  uncertainty  is  to 
achieve  a  90$  confidence  interval  with  a  halfwidth  of  0.1 
centered  on  the  P^.  estimate.  It  follows  then  that  an  error 
of  this  magnitude  resulting  from  the  use  of  a  characterization 
is  also  acceptable. 

In  terms  of  useful  results,  the  actual  characterization 
arrived  at  in  this  thesis  only  serves  to  demonstrate  the  feasi¬ 
bility  of  the  method.  The  behavior  of  any  particular  endgame 
program  is  dependent  upon  the  warhead  description  and  vulnera¬ 
bility  model  used  as  well  as  the  algorithm  employed  in  the 
program  itself.  What  should  be  of  interest  to  the  reader  is 
the  1116  thod  of  arriving  at  this  characterization,.  In  that  end¬ 
game  programs  are  all  models  of  the  same  basic  process,  this 
method  should  be  applicable  to  any  of  them.  Hence,  development 
of  the  method  of  obtaining  an  endgame  program  characterization 
is  the  most  important  objective  of  this  thesis. 


Utility  of  Results 


No  characterization  can  predict  as  well  as  the  endgame 
program  itself.  One  might  then  logically  question  the  wisdom 
of  pursuing  such  an  endeavor.  The  fact  is  that  a  characteri¬ 
zation  which  meets  the  stated  tolerance  criteria  for  a  parti¬ 
cular  endgame  program  would  have  several  potential  benefits. 

Once  established  and  verified,  the  characterization 
could  be  used  in  place  of  the  endgame  program/vulnerability 
model/warhead  description  combination  on  which  it  is  based. 

It  should  be  able  to  estimate  an  individual  Pj,  in  a  fraction 
of  a  central  processor  second.  In  addition,  the  characteri¬ 
zation  should  contribute  to  better  understanding  of  the 
processes  and  interrelationships  which  affect  the  endgame 
s i tuat iuh . 

There  are  several  applications  of  endgame  analysis 
which  would  also  benefit  from  a  characterization.  One  of 
these  is  real  time  missile  simulation  such  as  that  used  on 
the  Air  Combat  Maneuvering  Instrumentation/Range  (ACMI/R)  at 
Nellis  Air  Force  Base,  Nevada.  Currently,  the  results  of  the 
missile  simulation  are  used  to  enter  a  table  from  which  a  Pj. 
value  is  obtained.  The  increments  on  the  table  parameters 
are  relatively  large  and,  hence,  in  general  yield  a  crude  esti 
mate  of  P^.  Since  endgame  programs  are  used  to  generate  thes 
tables,  a  cliaracteri zation  based  on  the  same  amount  (or  less) 
of  data  as  is  in  the  tables  would  be  a  more  accurate  represen 
tation  of  the  program's  behavior,  would  use  less  computer 
memory  than  large  look-up  tables,  and  could  conceivably  be  a 
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faster  way  to  get  a  Pj,  estimate. 

A  related  application  is  in  the  area  of  missile  launch 
envelopes.  These  envelopes  are  used  by  pilots  to  help  them 
fire  their  missiles  under  conditions  which  supposedly  will 
result  in  a  "kill."  Currently,  these  envelopes  are  based  on 
miss  distance  or  some  other  criteria.  The  concept  of  using 
a  Pj,  criterion  has  long  been  proposed,  but  not  pursued  due  to 
the  inordinate  cost  of  the  endgame  analysis  required.  With 
a  characterization  of  the  endgame  program's  behavior,  such 
Pj.  envelopes  would  be  much  easier  and  cheaper  to  produce  and 
could  conceivably  even  be  implemented  in  aircraft  fire  control 
computers . 

A  final  application  of  endgame  analysis  which  could 
benefit  from  characterization  is  the  area  of  fuzing  analysis. 
Of  primary  importance  to  the  fuze  is  the  question  of  how  long 
the  time  delay  between  target  detection  and  warhead  detonation 
should  be.  Since  a  characterization  could  very  quickly  yield 
the  location  of  maximum  P^  for  any  given  trajectory,  one  need 
only  relate  this  location  to  the  fuze  detection  location  for 
the  same  trajectory.  Currently,  most  endgame  programs  must 
"search"  for  the  point  of  maximum  P^.  With  "searching"  being 
expensive  and  time-consuming,  such  analysis  has  been  limited 
in  the  past. 

Approach 

The  approach  to  the  problem  taken  herein  is  a  heuristic 
one.  Of  utmost  importance  is  the  availability  of  an  endgame 
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program  which  is  typical  of  those  in  use.  The  Missiles  and 
Guns  Analysis  Branch  of  the  Air  Force  Armament  Laboratory  has 
provided  a  version  of  its  primary  endgame  program  known  as 
SILAZAM.  The  warhead  description  and  vulnerability  model 
included  in  this  version  are  representative  of  the  type 
encountered  in  modern  endgame  analyses.  Hence,  the  program 
should  provide  a  realistic  challenge  in  terms  of  behavior 
characterization. 

Given  the  endgame  program,  the  approach  is  to  use  any 
available  methodology  to  establish  and  characterize  its  beha¬ 
vior.  To  the  extent  that  available  methodology  is  inadequate, 
it  is  supplemented  with  appropriate  "ad  hoc"  methods  to  accom¬ 
plish  the  same  end. 

Overview 

The  remainder  of  this  thesis  is  devoted  to  the  pursuit 
of  the  objectives.  The  second  chapter  contains  a  discussion  of 
available  methodologies  which  are  of  potential  usefulness  in 
this  problem.  Chapter  III  contains  a  description  of  the 
characterization  effort,  while  Chapter  IV  delineates  the 
resulting  suggested  procedure  for  obtaining  a  character! zat ion . 
The  final  chapter  contains  a  summary,  conclusions,  and  recom¬ 
mendations  for  further  research. 
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1 1 .  Available  Methodologies 


The  basic  problem  addressed  in  this  thesis  bears 
resemblance  to  several  which  have  been  treated  in  the  litera¬ 
ture.  In  fact,  there  are  three  somewhat  overlapping  "disci¬ 
plines"  which  are  concerned  with  the  general  problem  of 
describing  the  relationship  between  several  independent 
variables  and  a  dependent  variable  when  the  true  relationship 
is  either  not  known,  not  understood,  or  not  practically  useful. 

The  first  and,  perhaps,  oldest  of  these  disciplines  is 
numerical  analysis.  A  fundamental  concept  in  numerical 
analysis  is  the  collocation  polynomial.  It  can  be  mathemati¬ 
cally  shown  that  for  any  collection  of  observations  taken  at 
different  values  of  the  independent  variables,  there  exists  a 
unique  polynomial  which  takes  on  the  corresponding  values  of 
(i.e.,  collocates  with)  the  value  of  the  dependent  variable 
(Ref  2:10).  Much  of  numerical  analysis  is  devoted  to  methods 
for  finding  or  approximating  the  collocation  polynomial . 

The  main  disadvantage  of  collocation  polynomials  is 
that  there  is  no  assurance  regarding  the  behavior  (or  mis¬ 
behavior)  of  the  function  between  the  observed  points. 
Furthermore,  because  each  set  of  data  has  a  unique  polynomial 
which  collocates  with  it,  it  is  indeed  unlikely  that  two  inde¬ 
pendent  sets  of  observations  of  the  same  process  will  yield 
the  same  collocation  polynomial.  Finally,  the  order  of  the 
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collocation  polynomial  increases  with  the  number  of  observa¬ 
tions.  making  it  unwieldy  for  large  data  sets. 

Another  method  of  numerical  analysis  which  remedies 
some  of  these  problems  is  the  so-called  "spline  fit."  This 
involves  "piecing"  together  a  set  of  cubic  polynomials,  each 
of  which  fits  the  data  in  some  sub- interval .  At  the  points 
where  the  curves  join,  their  first  and  second  derivatives  are 
forced  to  be  equal.  This  method  is  most  useful  for  interpola¬ 
tion  and  for  calculating  integrals  and  derivatives.  It  is, 
however,  somewhat  cumbersome,  particularly  for  large  volumes 
of  data  where  more  than  one  independent  variable  is  involved 
(Ref  3  :474-488)  . 

Linear  regression  analysis  is  an  approach  more  suited 
to  large  volumes  of  data.  Rather  than  requiring  collocation, 
regression  analysis  seeks  to  minimize  the  error  between  the 
regression  function  and  the  observed  data.  The  regression 
function  is  a  linear  combination  of  different  functions  and 
must  be  linear  in  all  coefficients.  Transformations  can  often 
be  performed  on  both  the  data  and  the  proposed  function  in  order 
to  meet  these  requirements.  Given  the  data  and  the  proposed 
function,  linear  regression  analysis  yields  values  for  the 
coefficients  which  minimize  the  sum  of  the  squared  differences 
between  the  observed  values  of  the  dependent  variable  and  the 
corresponding  values  of  the  regression  function.  The  resulting 
equation  is  guaranteed  to  give  the  best  fit  (in  the  least 
squares  sense)  of  all  possible  equations  with  the  form  of  the 
regression  function.  There  is  no  guarantee,  however,  that 
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this  form  is  the  best  from  which  to  choose  (Ref  4;  5;  and  6). 

The  newest  discipline  which  deals  with  the  problem  is 
Response  Surface  Methodology  (RSM) .  On  the  surface,  at  least, 
it  appears  to  be  a  panacea  for  just  the  type  of  problem  repre¬ 
sented  by  this  thesis.  In  the  context  of  RSM,  one  has  a  process 
with  a  single  measurable  output  (response)  and  perhaps  several 
(say  m)  inputs  (factors) .  For  any  combination  of  values  of 
the  factors,  there  is  but  a  single  response.  Though  the  measure¬ 
ment  of  the  response  may  be  subject  to  error,  the  underlying 
process  can  be  viewed  as  a  surface  in  m  +  1  space.  RSM  seeks 
to  explore  this  surface  in  such  a  way  as  to  accomplish  some 
goal.  There  are  two  goals  with  which  RSM  is  concerned:  the 
attainment  of  optimum  conditions  and  the  characterization  of 
the  surface  itself.  Unfortunately,  the  overwhelming  majority 
of  the  literature  is  devoted  to  the  former  goal,  while  the 
treatment  of  the  latter  is  largely  lip-service. 

A  large  part  of  the  methodology  of  RSM  is  actually 
numerical  analysis  and  regression  analysis.  However,  a  signi¬ 
ficant  portion  of  the  discipline  is  also  devoted  to  experimental 
design- - i . e . ,  how  best  to  "examine"  the  surface  to  accomplish 
the  goal  with  a  minimum  of  experimentation.  This  portion  of 
the  discipline  is  of  particular  interest  to  anyone  faced  with 
the  basic  RSM  problem  (Ref  7;  8;  and  9). 

The  purpose  of  this  rather  cursory  discussion  of 
available  methods  has  been  to  acquaint  the  reader  with  them 
and  to  acknowledge  their  existence.  They  are  kept  in  the 
"characterization  tool  box"  for  use  as  needed. 
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III.  Analysis 


Factor  Selection 

In  any  endeavor  in  which  one  is  exploring  the  relation¬ 
ship  between  an  output  and  numerous  inputs ,  the  selection  of 
the  set  of  inputs  to  be  studied  is  of  utmost  importance.  Most 
processes  can  be  described  by  several  different  sets  of  factors 
(or  inputs) .  These  sets  are  generally  related  to  each  other 
somehow  and  might  be  overlapping,  but  they  do  not  necessarily 
contain  the  same  number  of  factors.  The  goal  is  to  select 
that  set  of  factors  which  is  most  suitable  for  use  in  predicting 
the  response  (or  output).  In  general,  the  factors  in  this  set 
have  the  most  direct  impact  upon  the  response  and  the  least 
amount  of  interaction  between  them.  If  the  researcher  is 
completely  ignorant  of  the  process  under  study,  a  large  amount 
of  experimentation  may  be  required  to  find  the  most  appropriate 
factor  set.  On  the  other  hand,  he  may  have  sufficient  knowledge 
of  the  process  to  choose  the  factor  set  without  any  experimen¬ 
tation  . 

In  the  case  of  the  endgame  situation,  there  are  several 
sets  of  variables  which  can  precisely  define  it.  There  are 
fundamentally  only  two  factors  involved:  spatial  relationship 
and  motion.  In  three  dimensions,  three  is  the  minimum  number 
of  parameters  required  to  describe  translational  motion.  The 
position  of  one  point  with  respect  to  another  point  in  three 
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dimensions  also  requires  a  minimum  of  three  descriptive  para¬ 
meters.  Additional  parameters  are  required  to  describe  the 
angular  orientation  and  motion  of  any  rigid  bodies  involved. 

Figure  2  depicts  that  part  of  the  geometry  which  results 
from  the  motion  of  the  target  and  the  missile.  The  coordinate 
system  is  the  target  coordinate  system  with  the  X-axis  coinci¬ 
dent  with  the  target  centerline,  the  Y-axis  directed  out  the 
port  side  of  the  target,  and  the  Z-axis  directed  out  of  the 
top  of  the  target  perpendicular  to  the  wing  plane.  The  origin 
of  the  system  is  the  target  centroid- - i . e .  the  geometric  center 
of  the  target.  For  purposes  of  this  study,  it  is  assumed  that 
the  target  centerline  is  aligned  with  the  target  velocity  vector 
and  that  the  target  is  not  rolled  out  of  the  horizontal  plane. 

In  the  definitions  presented  in  Table  1,  azimuth  refers  to  the 
angle  as  measured  in_  the  horizontal  (X-Y)  plane  and  elevation 
refers  to  the  angle  as  measured  from  the  horizontal  plane.  Also 
to  maintain  mathematical  correctness  ,  the  angles  in  these 
definitions  must  be  described  with  respect  to  the  negative  of 
both  the  missile  and  relative  velocity  vectors.  This  can  be 
reconciled  with  Figure  2  by  remembering  that  all  vectors,  by 
necessity,  are  defined  with  respect  to  the  origin  of  the 
coordinate  system.  In  Figure  2,  vectors  have  been  moved  from 
the  origin  in  order  to  more  clearly  illustrate  the  geometry. 

The  reader  should  refer  to  Figure  2  when  reading  the  definitions 

Figure  3  depicts  that  part  of  the  geometry  resulting 
from  the  spatial  relationships  between  the  missile  and  the 
target.  The  definitions  in  Table  11  apply  to  this  figure. 


Kinematic  Geometry 


TABLE  I 

Definitions  Related  to  the 
Kinematic  Geometry 


Term 

Definition 

the  target  velocity  vector 

\ 

the  missile  velocity  vector 

?R 

the  relative  velocity  vector  defined  by 
the  vector  difference  between  Vx.  and  V_ 
i  e  V  =  7  -  V  M  T 

x.e.  vR  vT 

VT 

magnitude  of 

VM 

magnitude  of  Vjj 

VR 

magnitude  of  VR 

azimuth  angle  of  -V..  with  respect  to  V_ 

•y 

elevation  angle  of  -V^  with  respect  to  VT 

'I' 

azimuth  angle  of  -VR  with  respect  to  VT 

X 

elevation  angle  of  with  respect  to  VT 

a 

angle  between  and  the  missile  centerline 

B 

angle  between  7R  and  the  missile  centerline 

t 

l 


I  l 

Fig  3.  Spatial 
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TABLE  II 


Definitions  Related  to  the 
Spatial  Geometry 


Term 

Definition 

Relative 

Traj  ectory 

path  of  the  missile  with  res¬ 
pect  to  the  target.  Direction 
corresponds  to  the  direction 
of  VR. 

Point  of 

Closest  Approach 

Point  on  the  relative  trajec¬ 
tory  at  which  the  line-of- 
sight  to  the  target  centroid 
is  perpendicular  to  VD 

Plane  of 

Closest  Approach 

Plane  which  contains  the  target 
centroid  and  is  perpendicular 
to  VR 

Relative 

System 

Coordinate 

Coordinate  system  obtained  by 
rotating  the  target  coordinate 
system  first  through  the  sup¬ 
plement  of  4;  and  then  through 
-X  such  that  the  X-axis  is 
aligned  with  VR.  Mote  that 
the  Y-Z  plane  is  the  plane 
of  closest  approach. 

XM(Burst 

point) 

X  coordinate  in  the  relative 
coordinate  system  of  the  point 
at  which  the  warhead  is 
detonated 

ym 

Y  coordinate  in  the  relative 
coordinate  system  of  the  point 
of  closest  approach 

ZM 

Z  coordinate  in  the  relative 
coordinate  system  of  the  point 
of  closest  approach 

R..(ri'iss 

M 

e 

distance) 

distance  between  point  of  clos¬ 
est  approach  and  target  centroid 

Angle  whose  tangent  is  YrraXA 
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Note  that  in  the  plane  of  closest  approach,  and  are 
the  rectangular  coordinates  of  the  point  of  closest  approach, 
while  R^  and  e  are  its  polar  coordinates. 

Thus,  to  define  a  particular  endgame  situation,  one  can 
specify  the  values  of  VM,  VT,  4C,  y,  a,  XM,  YM,  and  ZM,  or  one 
can  specify  the  values  of  VR,  4,  X,  6,  X^,  R^,  and  6,  or  one 
can  specify  some  other  set  of  parameters  which  contains  the 
essential  information.  To  limit  the  scope  of  this  project, 
the  value  of  y  (and  hence,  also  A)  is  assumed  to  be  zero. 
Therefore,  all  endgame  encounters  considered  herein  can  be 
completely  described  without  consideration  of  any  elevation 
angles.  For  purposes  of  methodology  demonstration,  this 
limitation  has  negligible  impact  because  the  means  of  account¬ 
ing  for  A  should  be  essentially  the  same  as  that  for  4 .  The 
assumption  is  made  merely  to  limit  the  size  of  the  project. 

The  set  of  factors  containing  VR,  g,  X^,  R^ ,  and  e 
is  deemed  the  most  appropriate  set  of  inputs  for  purposes  of 
studying  the  input/output  relationships  of  SHAZAM  for  several 
reasons.  First,  it  combines  the  information  provided  by  V^, 
VM,  4^ ,  and  a  into  just  three  factors:  VR ,  ^ ,  and  g. 

Secondly,  the  individual  effect  of  on  PR  depends  not  only 
on  the  value  of  ,  but  on  the  value  of  ^  as  well.  The  rela¬ 
tionship  of  PR  to  both  VR  and  4,  however,  is  more  direct  and 
less  dependent  on  interaction.  Finally,  the  effects  of  Y^ 
and  ZM  are  also  highly  interdependent,  while  the  use  of 
and  6  separates  the  effects  of  distance  and  orientation. 

Unfortunately,  SHAZAM  is  designed  to  accept  only  V^, 
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VT>  ipc ,  a,  Xjvj,  ,  and  2^  as  inputs  to  describe  the  geometry. 
The  relationship  between  these  inputs  and  those  mentioned 
above  is  given  in  Appendix  A. 

The  Empirical  Approach 

In  a  purely  empirical  approach,  it  is  necessary  to 
design  an  experiment  which  exposes  the  effects  of  the  factors 
on  Pj,,  There  are  both  individual  effects  and  combined  effects 
v/hich  must  be  accounted  for.  Classical  designs  such  as  those 
utilized  in  RSM  generally  require  a  large  number  of  data 
points.  For  example,  a  full  factorial  experiment  with  six 
factors  at  four  levels  each  requires  4096  data  points.  The 
amount  of  computer  time  required  by  SHAZAM  to  compute  a  single 
Pj,  is  dependent  upon  the  number  of  iterations  required  to 
achieve  the  specified  variance.  The  smallest  acceptable 
number  of  iterations  is  30,  while  no  more  than  70  iterations 
are  required  to  attain  a  90  percent  confidence  interval  about 
the  Pj,  estimate  with  a  halfwidth  of  0.1.  Under  these  condi¬ 
tions,  SHAZAM  requires  an  average  of  about  10  CP  seconds  (on 
a  CYBER  175)  to  compute  a  single  Pj,.  Hence,  a  full  factorial 
experiment  would  require  40,960  seconds  (or  more  than  11  hours) 
of  CP  time ! 

A  second  problem  with  classical  designs  is  that  in  the 
case  of  the  endgame  situation,  not  all  combinations  of  factor 
levels  would  yield  a  viable  endgame  encounter.  For  example, 
if  \p  -  180°,  B  must  be  very  small  in  order  for  the  geometry  to 
resemble  reality.  However,  for  =  90°,  3  cannot  be  very  small 
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at  all.  Hence,  many  data  points  would  have  to  be  omitted 
from  even  a  partial  factorial  experiment.  The  result  would 
be  the  loss  of  most  of  the  advantages  that  these  designs  are 
supposed  to  provide. 

Finally,  X^  exhibits  properties  of  both  a  factor  and 
a  response.  For  any  combination  of  the  other  five  factors, 
changing  X^  (the  burst  point)  can  cause  Pj,  no  change  from  0 
to  1.0  or  vice  versa.  Furthermore,  the  values  of  X..  which 
yield  non- zero  values  of  Pj,  are  highly  dependent  upon  the 
settings  of  the  other  factors. 

In  view  of  these  difficulties,  it  is  necessary  to  cus¬ 
tom  design  an  experiment  to  yield  the  necessary  information 
with  a  reasonable  amount  of  computer  time.  The  selected 
region  of  experimentation  is  presented  in  Table  III.  The 
range  of  values  for  VR  is  based  on  values  which  are  typically 
encountered  in  most  missile  systems.  Because  the  target  is 
symmetric  and  X  =  0 °,  need  only  be  examined  from  0°  (head-on) 
to  180°  (tail-on) .  The  value  of  6  can  range  from  0  to  an 
upper  limit  which  is  determined  by  the  maximum  look-angle  of 
the  missile's  seeker.  The  value  of  45°  is  fairiy  typical  of 
this  limit  for  many  missiles.  The  range  selected  for  RM 
encompasses  all  values  of  conceivable  interest  since  the  P^ 
at  100  feet  is  generally  very  close  to  zero  under  any  condi¬ 
tions.  Of  course,  9  can  have  any  value  from  0°  to  360°,  and 
the  value  of  X^  is  dependent  upon  the  other  factors  as  pre¬ 
viously  mentioned. 


TABLE  III 


Region  of  Experimentation 


Factor 

Minimum 

Maximui, 

VR 

1000  ft/sec 

3000  ft/sec 

0° 

180° 

B 

0° 

45c 

*M 

0  ft 

100  ft 

e 

0° 

360° 

XM 

as  needed 

as  needed 

Table  IV  is  a  staged  experimental  design  over  this 
region.  The  underlying  premise  is  that  the  highest  order  of 
significant  factor  interaction  is  three.  The  first  stage  is 
designed  to  yield  information  about  the  individual  and  com¬ 
bined  effects  of  tjj ,  e,  and  X^.  The  second  stage  is  designed 
to  expose  any  change  in  these  effects  if  only  one  of  the  re¬ 
maining  three  factors  is  changed.  Stage  three  examines  the 
effect  of  simultaneously  changing  any  two  of  the  factors  VR, 

,  and  g.  Note  that  in  stage  three,  only  one  combination 
of  ip  and  0  is  called  for.  The  assumption  is  that  the  effects 
discovered  in  this  stage  are  the  same  for  all  combinations  of 
ijj  and  e.  In  stage  four,  the  effect  of  simultaneously  varying 
all  three  of  VR,  R^ ,  and  g  is  examined.  Again,  this  effect 
is  assumed  to  be  independent  of  ip  and  6.  Stage  five  is  simply 
a  means  for  verifying  that  the  assumptions  hold. 

For  reasons  which  are  explained  later,  only  the  first 
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TABLE  iV 


Staged  Experimental  Design 


r 

e° 

V£t) 

VR(fps) 

vft) 

Stage 

1  : 

0 

45 

45 

135 

90 

225 

5 

2000 

15 

15 

135 

315 

NZYV* 

180 

! 

Stage 

2: 

45 

45 

5 

1000 

15 

15 

180 

225 

NZYV 

3000 

4  5 

45 

5 

2000 

27.5 

15 

180 

225 

NZYV 

40 

45 

45 

5 

2000 

15 

7  .  5 

180 

225 

NZYV 

22.5 

Stage 

_3 : 

135 

225 

5 

1000 

27 , 5 

15 

NZYV 

3000 

40 

135 

225 

5 

1000 

15 

7.5 

NZYV 

3000 

22.5 

135 

225 

5 

2000 

27  .  5 

7.5 

NZYV 

40 

22.5 

Stage 

A'- 

135 

225 

5 

1000 

27  .  5 

7.5 

NZYV 

3000 

40 

22.5 

Stage 

_5  : 

Experiment  as 

needed. 

*NZYV 

=  Non  - 

zero  yielding 

values 

stage  of  this  experiment  was  completed.  The  following  is  a 
discussion  of  the  analysis  of  the  data  generated  by  this  first 
stage . 

For  each  combination  of  ^  and  e,  the  data  were  plotted 

with  P,,  as  the  ordinate  and  X.,  as  the  abscissa.  All  of  these 
X  M 

plots  displayed  a  common  characteristic:  in  each  case,  P„ 
was  zero  until  a  very  steep  rise  to  one,  where  it  stayed  for 
some  distance  and  then  rapidly  dropped  to  zero  again.  The 
width  and  location  of  this  "pulse"  varied  with  \p  and  0. 

All  attempts  to  fit  polynomial  functions  to  these  pulses 
using  either  numerical  analysis  or  linear  regression  were  fruit¬ 
less.  No  polynomial  of  reasonably  low  order  could  be  found  to 
fit  all  of  the  pulses.  In  fact,  to  fit  any  single  pulse,  it 
was  necessary  to  "clip"  every  polynomial  that  was  examined  to 
prevent  it  from  either  exceeding  one  or  going  below  zero.  It 
became  apparent  that  some  other  means  of  fitting  these  pulses 
was  needed. 

After  numerous  candidate  functions  were  examined,  it 
was  found  that  the  function: 

PK  -  e'(-V) 

could  be  made  to  fit  each  plot  reasonably  well  just  by  changing 
a  and  b.  In  addition,  no  clipping  was  necessary  since  the 
function  can  never  exceed  one  nor  be  less  than  zero.  Unfortun¬ 
ately,  the  non-linear  nature  of  this  function  virtually  elimin¬ 
ates  RSM  as  a  useful  tool. 
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The  problem  then  became  one  of  accounting  for  the  changes 
in  a  and  b  due  to  changes  in  \p  and  9.  The  value  of  a  corres¬ 
ponds  to  the  negative  of  the  value  of  X^  at  the  center  of  the 
pulse  while  b  is  the  halfwidth  of  the  pulse  as  measured  between 
the  points  =  1/e  on  the  up  and  down  side. 

For  each  data  plot,  the  values  of  a  and  b  were  visually 
estimated.  These  estimates  were  plotted  versus  0 ,  holding 
constant.  In  each  case,  it  was  found  that  a  function  of  the 
form 

c1  *  c2  sin2  (|e  f  C3) 


could  be  used  to  reasonably  fit  each  of  these  plots  (see 
Appendix  B,  Figures  1  through  5  for  a  and  Figures  6  through  10 
for  b)  .  It  then  remained  to  describe  the  changes  in  C^,  , 

and  C3  i  ith  changes  in  tj;  for  the  "location"  (a)  function  and 
the  "halfwidth"  (b)  function  (see  Appendix  B,  Figures  11  through 
13  for  location  and  Figures  14  through  16  for  halfwidth) .  The 
resulting  functions  were  as  follows: 


\  /i  / 


?  Q 


a  =  5.5  +  5  sin<J>  -  10  \f  sin(jgg-)  sin"(^0 


10 


5  sin  (^jj)  -  5e 


,-\p  -  100^8  ru-. 

"CT4 - }  [Hj 


where 


h 


;in2Ly  6  +  90cos2  (1---3^(^-— --}-)  cos3 ( 1 . 33 -  45 j ) ] 


i 

i 

t 


Thus,  the  function  arrived  at  to  predict  P^,  given  if;, 
0 ,  and  X^  was : 
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with  a  and  b  as  described  above.  For  each  of  the  16  combina¬ 
tions  of  4)  and  6  generated  in  stage  one,  this  function  was 
plotted  along  with  the  data  points  to  visually  examine  the  fit 
(see  Appendix  B,  Figures  17  through  32). 

Examination  of  these  plots  led  to  several  conclusions. 
First,  although  the  basic  behavior  of  SHAZAM  was  captured  by 
the  function,  there  was  substantial  error  on  some  of  the  plots. 
This  error  was  partially  due  to  the  fact  that  in  some  cases, 
the  rise  was  either  too  steep  or  too  shallow.  A  second  con¬ 
tributor  to  the  error  was  the  generally  poor  fits  to  the 
functions  for  C-^,  C.^,  C„  ,  a,  and  b.  Obviously  a  second  itera¬ 
tion  was  required  using  linear  regression  to  fit  these  functions 
rather  than  rough  calculation. 

Ihe  second  conclusion  was  that  while  the  functions  for 
a  and  b  only  accounted  for  two  of  the  factors  4,  9,  V^,  R^ , 
and  8,  they  had  already  become  immensely  cumbersome  and  almost 
intractable.  How  much  more  so  would  they  become  for  all  five 
factors? 

In  contemplation  of  this  question,  the  conclusion  was 
reached  that  the  original  assumption  concerning  the  insignifi¬ 
cance  of  higher  than  third  order  interactions  was  incorrect. 
Clearly,  since  4,  e,  V^,  ,  and  B  ate  all  essential  in  des¬ 

cribing  the  endgame  geometry,  then  the  effect  of  any  one  of 
these  factors  on  the  width  and  location  of  the  pulse  is 
dependent,  on  the  settings  of  all  four  of  the  remaining  factors. 
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To  discover  the  functional  relationship  which  accoi  its  for 
this  fifth  order  interaction  would  require-  an  excessive 
amount  of  experimentation.  Whether  such  a  function  could  be 
found  is  questionable,  and  even  if  it  could,  it  would  prob¬ 
ably  be  so  unmanageable  and  error  prone  that  it  would  be  of 
limited  usefulness.  Furthermore,  in  terms  of  a  general 
methodology,  it  must  be  remembered  that  an  additional  factor, 
elevation  angle,  must  eventually  be  included.  Therefore,  it 
was  decided  to  abandon  the  purely  empirical  approach  in  favor 
of  a  method  utilising  the  geometry  itself  to  account  for  the 
interactive  effects  of  the  factors  on  a  and  b. 

A  Semi -Empirical  Approach 

The  idea  in  this  approach  is  to  explain  as  much  as 
possible  of  the  variation  in  a  and  b  caused  by  the  factors 
,  VR,  Rm,  6,  and  B  with  the  geometry  itself.  The  remaining 
variation  is  then  explained  on  the  basis  of  empirical  obser¬ 
vation.  Hence,  this  semi-empirical  approach  is  a  combination 
of  analytical  and  empirical  methods. 

Recall  that  the  purpose  of  an  endgame  program  is  to 
account  for  the  interaction  between  the  warhead  characteristics, 
the  geometry,  and  the  vulnerability  model.  This  interaction 
can  be  viewed  as  consisting  of  two  processes:  the  interactive 
effect  of  the  warhead  characteristics  and  the  geometry  on  the 
motion  of  the  warhead  fragments,  and  the  interaction  between 
the  fragments  and  the  vulnerability  model.  Note  that  blast 
effects  have  been  excluded  in  this  discussion  and  will  be 
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addressed  later. 

The  interaction  between  the  fragments  and  the  vulnera¬ 
bility  model  is  determined  by  the  nature  of  the  fragments 
themselves  as  well  as  by  the  direction  and  magnitude  of  their 
velocity  vectors.  None  of  the  parameters  which  describe  the 
geometry  have  any  effect  upon  such  things  as  fragment  mass 
and  size.  ,  However,  the  geometry  does  have  a  profound  effect 

4 

on  the ‘speed  and  direction  of  travel  of  the  fragments. 

The  interaction  between  the  warhead  characteristics 
and  the  geometry  can  be  analytically  calculated  given  enough 
information  about  the  warhead  model.  The  information  needed 
concerns  the  direction  and  speed  with  which  fragments  are 
ejected  from  the  warhead  under  static  conditions.  For  the 
warhead  model  in  the  version  of  SHAZAM  used  in  this  study, 
the  fragment  velocity,  Vp ,  is  5500  feet  per  second.  There 
are  two  directions  between  which  fragments  are  ejected:  84° 
and  91°  as  measured  from  the  missile  centerline. 

An  accurate  way  to  describe  fragment  motion  is  to  con¬ 
sider  each  fragment  individually.  In  this  method,  the  fragment 
velocity  vector  is  added  to  the  relative  velocity  vector  to 
determine  the  path  of  the  fragment  relative  to  the  target. 

Another  way  of  accounting  for  fragment  motion  is  to 
consider  the  static  fragment  pattern  as  a  cone  emanating  from 
the  warhead.  In  the  dynamic  situation,  the  orientation  of 
the  cone  with  respect  to  tie  missile  centerline  as  well  as 
the  half-angle  of  the  cone  are  changed  according  to  the  magni¬ 
tude  and  direction  of  the  relative  velocity  vector.  This 
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process  is  an  approximation  to  the  true  dynamic  warhead  pattern 
described  in  the  preceding  paragraph.  The  approximation  is 
generally  very  good,  but  some  error  is  introduced  in  the  case 
when  3  is  large  and  the  conic  section  described  by  the  inter¬ 
section  of  the  X-Z  plane  (in  the  target  coordinate  system)  and 
the  cone  is  an  ellipse  or  a  circle.  The  method  of  calculating 
the  cone  and  some  illustrations  of  the  approximation  are  given 
in  Appendix  C . 

The  advantage  of  a  cone  description  of  the  warhead  pattern 
is  that  the  intersection  between  the  cone  and  the  infinite  line 
containing  the  target  centerline  can  be  obtained  analytically. 
There  are  two  ways  to  solve  for  this  intersection:  1)  given 
a  specific  point  on  the  line,  determine  the  value  of  X^  at 
which  the  point  lies  on  the  cone;  or,  2)  given  X^ ,  find  the 
point  or  points  on  the  line,  if  any,  which  are  also  on  the  cone. 
The  mathematical  solutions  to  these  two  problems  are  given  in 
Appendix  D. 

For  each  data  point  generated  in  stage  one  of  the  experi¬ 
ment  described  in  the  preceding  section,  the  points  on  the  tar¬ 
get  centerline  lying  on  the  cone  were  obtained  using  the  second 
method  of  solution  mentioned  above.  There  seemed  to  be  a 
relationship  between  where  the  pulse  "jumped"  and  where  the  cone 
hit  on  the  target  centerline.  In  fact,  it  seemed  that  whenever 
at  least  one  of  the  two  fragment  cones  (one  cone  representing 
the  leading  edge  and  one  the  trailing  edge  of  the  warhead  pattern) 
was  hitting  the  target  centerline  between  about  seven  feet  behind 
the  centroid  and  about  five  feet  in  front  of  the  centroid,  the 


Pj,  was  near  its  highest  value. 

To  further  test  this  relationship,  more  data  were 
generated  using  differing  values  of  and  R^.  In  all,  a 
total  of  48  Pj,  versus  data  plots  were  produced.  The  com¬ 
binations  of  factor  values  examined  are  given  in  Table  V. 
Examination  of  these  plots  revealed  that  the  relationship 
was  fairly  good  except  for  situations  when  was  large.  In 
such  cases,  the  maximum  value  of  seemed  to  decrease,  while 
the  "vulnerable  length"  along  the  target  centerline  seemed  to 
shrink  considerably.  While  the  decrease  in  P^  can  be  ration¬ 
alized,  the  shrinkage  in  the  vulnerable  length  seemed  excessive. 
The  pulse  became  exceedingly  narrow  for  only  moderately  high 
values  of  R, ^ . 

A  possible  reason  for  this  apparent  shrinkage  in  the 
vulnerable  length  was  the  use  of  the  random  angle -of - attack 
(«)  option  in  SHAZAM  in  the  generation  of  the  data.  This  option 
was  intentionally  exercised  for  the  original  experiment  in  hopes 
of  making  the  rise  from  P^  -  0  to  maximum  Pj,  more  shallow  and, 
hence,  possibly  easier  to  fit.  In  fact,  for  larger  miss  dis¬ 
tances,  it  had  the  opposite  effect. 

To  understand  the  reason  for  this,  one  must  understand 
the  way  that  SHAZAM  treats  “  as  a  random  variable  when  that 
option  is  exercised.  A  mean  value  and  standard  deviation  for 
«  must  be  specified.  Then,  for  each  iteration,  SHAZAM  selects 
a  random  sample  from  a  normally  distributed  population  with 
mean  and  standard  deviation  as  specified.  This  value  is  then 
used  as  the  value  for  a  for  that  s ingle  iteration.  For  the 
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TABLE  V 


Combination  of  Factor  Values  Examined  in  Basic  Data 


6(°3 

VR 

(fps) 

Si 

(ft) 

B(°) 

<K°) 

e  (°) 

VR 

(fps) 

M 

(ft) 

B(°) 

0 

225 

2000 

15 

0 

135 

135 

2000 

27.5 

15 

0 

225 

2000 

27.5 

0 

135 

135 

2000 

40 

15 

0 

225 

2000 

40 

0 

135 

225 

1000 

15 

15 

0 

225 

3000 

15 

0 

135 

225 

2000 

15 

15 

0 

13  5 

2000 

15 

0 

135 

315 

2000 

15 

15 

0 

13  5 

3000 

15 

0 

180 

135 

100  0 

15 

0 

45 

45 

2000 

15 

15 

180 

135 

1000 

27.5 

0 

4  5 

45 

2000 

27.5 

15 

180 

135 

1000 

40 

0 

45 

4  5 

2000 

40 

15 

180 

135 

2000 

15 

0 

45 

135 

2000 

15 

15 

180 

135 

2000 

27.5 

0 

45 

135 

3000 

15 

15 

180 

135 

2000 

40 

0 

45 

225 

2000 

15 

15 

180 

135 

3000 

15 

0 

45 

315 

2000 

15 

15 

180 

135 

3000 

27 .  S 

0 

45 

315 

3000 

15 

15 

180 

135 

3000 

40 

0 

90 

45 

2000 

15 

15 

180 

225 

1000 

15 

0 

90 

135 

2000 

15 

15 

180 

225 

1000 

27.5 

0 

90 

135 

3000 

15 

15 

180 

225 

1000 

40 

0 

90 

225 

2000 

15 

180 

225 

2000 

15 

0 

90 

225 

2000 

27.5 

15 

180 

225 

2000 

27.5 

0 

90 

225 

2000 

40 

15 

180 

225 

2000 

40 

0 

90 

315 

2000 

15 

15 

180 

225 

2000 

52.5 

0 

135 

45 

1000 

15 

15 

180 

225 

3000 

15 

0 

135 

45 

2000 

15 

15 

180 

225 

3000 

27.5 

0 

135 

135 

2000 

15 

15 

180 

225 

3000 

40 

0 

next  iteration,  a  new  value  of  a  is  chosen.  The  effect  is 
that  for  each  burst  point,  the  missile  centerline  "oscillates" 
back  and  forth,  changing  from  iteration  to  iteration. 

For  all  of  the  data  points  thus  far  generated,  the  mean 
value  of  a  was  set  at  3°  while  the  standard  deviation  was  set 
at  1°.  This  had  little  effect  at  small  values  of  R.,  since  3’ 
results  in  a  deflection  of  only  52  feet  per  1000  feet  of  dis¬ 
tance.  Hence,  for  =  15  feet,  the  fragment  cone  moved  less 
than  a  foot  along  the  target  centerline.  However,  at  = 

50  feet,  the  fragment  cone  moved  up  to  2.6  feet!  This  dis¬ 
tance  was  sufficient  to  pull  the  fragment  cone  out  of  the 
vulnerable  length  completely.  The  result  was  that,  for  a 
burst  point  at  which  a  constant  a  =  0° would  have  yielded  a 
high  P~  due  to  the  leading  fragment  cone  intersecting  the 
vulnerable  length,  the  oscillations  of  a  moved  the  cone  away 
from  the  vulnerable  length  a  sufficient  number  of  times  to 
lower  the  P  .  A  comparison  between  data  points  generated  for 
the  same  trajectory  both  with  and  without  a  as  a  random  vari¬ 
able  is  shown  in  Figure  4.  Note  that,  as  expected,  setting 
a  =  0°  made  the  pulse  wider  . 

In  view  of  this  effect,  all  of  the  data  were  re-generated 
with  a  set  to  zero.  Some  shrinkage  was  still  observed  in  the 
vulnerable  length  as  became  large.  This  can  be  explained  by 
noting  that  since  the  fragment  pattern  expands  with  time,  the 
fragment  density  grows  smaller  and,  consequently,  the  pattern 
is  less  lethal.  This  explanation  can  also  be  applied  to  the 
decrease  in  maximum  Pj,  with  increasing  values  of  .  At  very 
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large  miss  distances  (>  75  feet) ,  the  vulnerable  length 
appeared  to  grow  again.  Since  miss  distances  of  these  magni¬ 
tudes  have  low  P^'s  and  are  generally  rare,  this  discrepancy 
has  little  impact  on  the  characterization. 

In  examining  the  newly  generated  data,  it  was  observed 
that  the  vulnerable  length  appeared  longer  when  6  was  between 
0°  and  180°  than  it  did  when  0  was  between  180°  and  360°.  The 
only  possible  explanation  for  this  is  that  the  target  is  appar¬ 
ently  more  vulnerable  from  the  top  than  it  is  from  the  bottom. 
This  also  manifests  itself  in  the  fact  that  maximum  tends 
to  decrease  more  slowly  with  increasing  for  trajectories 
above  the  target  as  opposed  to  those  below. 

The  importance  of  finding  this  "vulnerable  length" 
relationship  is  profound.  It  allows  a  and  b  to  be  determined 
strictly  from  the  geometry.  One  need  only  determine  what  pro¬ 
portion  of  maximum  P^  is  achieved  when  the  vulnerable  length 
is  first  contacted  by  a  fragment  cone.  Because  of  the  symmetry 
of  the  function,  the  same  proportion  of  maximum  Pj,  will  be 
achieved  when  the  vulnerable  length  is  last  contacted  by  a 
fragment  cone.  The  determination  of  a  and  b,  given  this 
information,  is  illustrated  in  Appendix  E. 

Unfortunately,  the  endpoints  of  the  vulnerable  length 
and  the  proportion  (P)  of  maximum  Pj.  achieved  at  these  points 
cannot  be  arrived  at  independently.  The  selection  of  each 
has  its  own  effect  on  the  resultant  pulse.  Therefore,  the 
selection  process  is  an  iterative  one- -different  combinations 
of  endpoints  and  proportions  must  be  tried  and  compared  to  the 
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data.  The  final  selection  of  these  values  is  strictly  a  matter 
of  the  experimenter's  judgment  as  to  which  combination  seems  to 
result  in  the  best  fit. 

In  pursuit  of  this  "final  selection,"  it  was  found  that 
a  fixed  vulnerable  length  would  not  suffice.  As  previously  men 
tioned,  the  vulnerable  length  tended  to  shrink  as  increased 
and  this  shrinkage  occurred  at  different  rates,  depending  on 
the  value  of  0.  Furthermore,  the  endpoints  of  the  vulnerable 
length  seemed  to  move  as  moved  from  0°  to  180°.  In  fact, 
for  a  given  value  of  P,  the  vulnerable  length  which  yielded 
the  best  fit  was  different  for  each  of  the  48  data  plots.  The 
problem,  then,  became  one  of  accounting  for  the  effects  of  the 
five  factors  on  the  endpoints  of  the  vulnerable  length. 

The  reason  for  the  changes  in  the  endpoints  of  the 
vulnerable  length  with  changes  in  the  factors  is  due  to  the 
structure  of  the  vulnerability  model.  Typically,  these  models 
are  very  detailed  descriptions  of  the  vulnerable  components  of 
the  target.  Fragment  damage  to  these  components  is  limited  in 
the  model  by  such  things  as  shielding,  angle  of  incidence, 
velocity,  mass,  and  fragment  density.  Obviously,  the  best  way 
to  capture  the  effects  of  all  of  these  parameters  is  with  the 
vulnerability  model  itself.  This,  however,  is  what  an  endgame 
program  does  and  is  obviously  the  cause  of  the  excessive  com¬ 
puter  time  requirement. 

For  the  approach  taken  herein,  the  vulnerable  length 
is  used  as  a  surrogate  for  the  vulnerability  model.  The 
ultimate  question  is:  how  much  of  the  behavior  of  the 
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vulnerability  model  must  be  accounted  for  with  changes  in  the 
vulnerable  length?  The  answer  to  this  question  is  again  a  matter 
of  the  experimenter's  subjective  judgment  as  to  what  constitutes 
an  adequate  fit.  Thus,  the  experimenter  must  observe  as  much 
of  the  behavior  as  possible  in  the  data  and  then  attempt  to 
account  for  that  which  he  or  she  thinks  is  important.  If  the 
results  are  then  unsatisfactory,  another  iteration  of  the  pro¬ 
cess  is  required  to  account  for  more  of  the  vulnerability 
model's  behavior. 

For  the  data  generated  in  this  project,  three  itera¬ 
tions  were  performed.  The  first  iteration  has  essentially 
already  been  discussed- -that  is,  the  discovery  of  a  vulnerable 
length  with  endpoints  seven  feet  before  and  five  feet  after 
the  target  centroid  that  seemed  to  provide  a  "rough'’  correla¬ 
tion  with  the  data.  In  the  second  iteration,  the  exponent 
in  the  function  was  changed  from  eight  to  six.  The  result 
of  this  was  to  make  the  pulse  rise  and  fall  more  slowly,  w'hich 
seemed  to  be  a  more  accurate  representation  of  the  data.  In 
addition,  a  linear  function  was  used  to  describe  the  decrease 
in  maximum  P^.  Also,  the  rear  endpoint  of  the  vulnerable  length 
was  moved  toward  the  forward  endpoint  as  a  linear  function  of 
R^.  The  details  of  these  adjustments  as  well  as  plots  of  the 
resultant  characterization  for  all  48  factor  combinations  are 
presented  in  Appendix  F.  The  plots  show  that  while  in  many  cir¬ 
cumstances  the  fit  is  very  good,  there  are  some  situations  where 
the  detail  in  the  vulnerability  model  was  not  captured  very  well 
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at  all.  The  latter  is  especially  true  for  the  larger  values 
of  Rjj. 

The  impact  of  these  errors  can  only  be  assessed  when 
the  results  are  compared  to  results  produced  by  SHAZAM  in  its 
normal  mode  of  operation.  Normally,  SHAZAM  (as  well  as  other 
endgame  programs)  treats  and  6  as  random  variables.  The 
reason  is  that,  in  reality,  the  miss  distance  achievable  by  a 
particular  missile  under  a  given  set  of  launch  conditions 
is  uncertain.  Flight  tests  over  the  years  have  shown  that 
although  other  parameters  of  the  endgame  geometry  can  be 
fairly  well  duplicated  from  firing  to  firing,  the  values  of 
R^  and  0  cannot  be  duplicated  with  any  certainty,  even  though 
the  firing  conditions  are  as  closely  duplicated  as  possible. 

The  reasons  for  this  phenomenon  are  too  numerous  and  complex 
to  discuss  here,  but  the  important  point  is  that  in  actual 
missile  endgame  analyses,  R^  and  0  (and,  consequently,  Y^  and 
Z^)  are  most  often  treated  as  random  variables.  The  most  com¬ 
mon  probability  distribution  used  to  describe  the  distribution 
of  miss  distance  is  the  bivariate  normal  (or  Gaussian)  distri¬ 
bution,  wherein  and  ZM  are  each  assumed  to  be  normally 
distributed  but  independently  of  each  other.  The  mean  of 
each  is  usually  taken  to  be  zero,  while  the  standard  deviations 
are  usually  made  equal  and  denoted  by  o.  The  result  is  what 
is  sometimes  referred  to  as  the  "Circular  Normal  Distribution." 
In  terms  of  R^  and  0 ,  the  equivalent  distribution  is  the  Ray¬ 
leigh  distribution.  An  important  point  about  this  distribution 
is  that  only  about  39  percent  of  the  possible  values  of  R^  lie 
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within  one  standard  deviation  from  the  mean,  whereas  roughly 
68  percent  lie  within  one  standard  deviation  from  the  mean 
in  the  univariate  normal  distribution. 

When  SHAZAM  treats  miss  distance  as  a  random  variable, 
values  for  Y^  and  2^  are  selected  from  a  normal  distribution 
in  a  psuedo-random  fashion.  A  new  selection  of  values  for 
and  ,  as  well  as  all  of  the  other  random  variables  in  the 
model,  is  made  for  each  Iteration.  As  previously  mentioned, 
the  total  number  of  samples  required  for  each  specification 
of  ip ,  V_,  6,  and  o  is  between  30  and  70  to  achieve  a  90  per¬ 
cent  confidence  interval  about  the  P.,  estimate  with  a  half- 

K 

width  of  0.1.  The  interpretation  of  this  confidence  interval 
is  that  there  is  a  .9  probability  that  the  estimate  is  within 
0.1  of  the  "true"  estimate.  By  "true"  estimate  is  meant  the 
estimate  that  would  result  if  the  number  of  samples  approached 
infinity . 

To  compare  the  results  of  the  characterization  to  those 
of  SHAZAM  under  these  circumstances,  it  is  necessary  to  account 
for  the  randomness  in  and  0..  What  is  actually  called  for 
is  the  following  integration: 

%„  -  Joro"  W»>  Vrh-U<Mrm  (3.1) 

where 

Pj  (Rji.e)  is  the  value  of  Fj.  at  and  0  and 

P  (Rm,0)  is  the  joint  probability  of  occurrence 
r  M  of  and  0 

With  the  characterization,  can  easily  be  found,  while 
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P;r  (^M ’  results  from  the  Rayleigh  distribution.  This  inte¬ 
gration  can  be  performed  numerically  with  what  is  known  as 
a  "109  Cell  Model”  of  the  Rayleigh  Distribution.  What  this 
model  does  is  to  brec  jp  the  plane  of  closest  approach  into 
109  cells.  For  97  of  the  cells,  the  probability  that  and 
6  lie  within  each  cell  is  .01.  For  the  remaining  12  cells, 
this  probability  is  .0025.  The  model  is  illustrated  in 
Appendix  G. 

One  w-'-y  of  estimating  Eq  (3.1)  using  the  characteriza¬ 
tion  is  to  find  values  of  a  and  b  for  each  of  the  109  relative 
trajectories  which  penetrate  the  probabilistic  "centers"  of 
each  of  the  cells.  Then 
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where  f .  (X,, , a  .  ,b  . )  is  the  characterizati  on  function  for  the 

1  M  x  1 

combination  of  factors  ip ,  V-> .  6,  R..  ,  and  9..  R..  and  0. 

i :  •  *  M .  l  M .  l 

i  th  1 

represent  the  probabilistic  "center"  of  the  i  cell. 

There  are  two  reasons  for  approaching  the  miss  distance 
distribution  problem  in  this  way.  First,  there  is  no  variance 
in  the  result.  The  procedure  will  produce  the  same  answer 
every  time.  The  erior  introduced  by  using  Eq  (3.2)  to  estimate 
F  ’  (3.1)  results  from  the  sum  of  the  weighted  differences 
between  the  Pj.  at  the  probabilistic  "center"  of  each  cell  and 
the  true  average  P^.  for  that  cell  as  predicted  by  the 


characterization.  The  weight  for  the  difference  is  just  the 
probability  associated  with  that  cell.  For  example,  if  the 
average  wi thin-cell  error  is  0.1,  then  the  total  error  induced 
by  the  109  Cell  Model  (E^)  is 

97  109 

E.,  =  I  (.011(0.1)  +  2  (.0025)  (0.1)  =  0.1  (3.3) 

M  i  =  l  i-9  8 

Since  most  of  the  cells  are  small,  one  would  expect  a  much 
smaller  average  within-cell  error.  Note  that  is  just  the 
error  caused  by  using  Eq  (3.2)  to  estimate  Eq  (3.1).  The 
error  in  using  the  characterization  instead  of  SHA2AM  to  esti¬ 
mate  P„  can  be  estimated  by 
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-  2  (.0025) (PK  -P„  )  (3.4) 

i  =  9  8  KCi  KSHi 


where  '  is  the  average  P^  given  by  the  characterization  for 
Ci 

each  cell  and  Pv  is  the  previously  mentioned  "true"  SHAZAM 

kSH  . 
i 

estimate  of  the  average  cell  P^. 

The  second  reason  to  use  the  "109  Cell  Model"  is  that 

for  a  given  combination  of  ty,  V^,  g,  and  one  need  only 

calculate  109  pairs  of  values  for  a  and  b.  Thereafter,  an 

infinite  number  of  values  of  X,,  can  be  examined. 

M 

Actual  calculation  of  using  Eq  (3.4)  is  impractical. 

A  total  of  109  values  of  P„  would  have  to  be  calculated. 

KSH 

Since  P,.  can  only  be  estimated  by  SHAZAM,  the  error  in  the 
KSH 

estimate  would  confound  the  true  value  of  each  of  the  (P». 

Cf 

P-.  )  terms.  Furthermore,  a  minimum  of  30x109  =  3270 

KSII. 
i 
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iterations  (1/2  hour  of  CP  time)  of  SHAZAM  would  be  required 
just  to  get  a  single  estimate  of  for  a  given  combination 
of  factor  values. 

Eq  (3.4)  does,  however,  give  some  indication  of  the 

sensitivity  of  E^  to  the  hinds  of  errors  observable  in  the 

plots  in  Appendix  F.  The  maximum  possible  error  in  any  cell 

is  1.0.  The  plots  show  that,  in  general,  the  difference 

between  the  SHAZAM  Pr  and  the  characterization's  is  much 

closer  to  (and  quite  often  less  than)  0.1.  Furthermore, 

since  some  of  the  Pv  's  are  less  than  the  corresponding  esti- 

KC 

mates  of  the  Pv  's  while  others  are  larger,  there  is  a 
KSH 

tendency  to  cancel.  Therefore,  it  seems  a  reasonable  expecta¬ 
tion  that  the  E^  of  the  characterization  is  on  the  order  of 
0.1  or  less.  If  Em  is  negligible,  then  a  total  error  of  0.1 
or  less  also  seems  a  reasonable  expectation.  Figure  5  is  an 
overlay  of  an  evaluation  of  Eq  (3.2) (for  a  candidate  charac¬ 
terization)  superimposed  on  an  estimate  of  Eq  (3.1)  based  on 
1000  samples  using  the  same  characterization.  It  appears 
that  E^  (the  difference  between  the  two  curves)  is  indeed 
very  small.  (The  "109  Cell"  curve  is  the  lower  one.) 

The  only  practical  way  to  estimate  E^  is  to  execute 
both  SHAZAM  and  the  characterization  for  identical  settings  of 
Rm,  ,  0,  and  0 ,  and  then  compare  the  results.  This  was  done 
for  13  such  combinations  using  the  characterization  shown  in 
Appendix  F.  The  lesulting  data  plots  are  shown  in  Appendix  H. 

It  is  appropriate  at  this  point  to  discuss  the  manner 
in  which  SHAZAM  deals  with  the  problem  of  direct  hit.  When 
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sampling  from  a  circular  normal  distribution,  there  is  a 
certain  probability  of  selecting  values  of  and  such 
that  the  relative  trajectory  penetrates  the  target  fuselage. 
The  version  of  SHAZAM  used  in  this  study  assigns  the  value  of 
one  to  the  for  all  such  trajectories  regardless  of  burst 
point.  In  order  to  match  the  SHAZAM  results,  it  is  necessary 
to  check  for  this  condition  before  evaluating  the  characteri¬ 
zation  function.  Without  knowing  any  of  the  physical  attri¬ 
butes  of  the  target,  it  was  necessary  to  experimentally 
determine  the  target's  length  and  width.  A  width  of  about  5 
feet  and  a  length  of  about  30  feet  seemed  to  match  the  SHAZAM 
direct  hit  calculations  reasonably  well.  Hence,  in  using  the 
109  Cell  Model,  the  trajectory  represented  by  each  cell  must 
first  be  checked  to  see  if  j |  is  less  than  2.3  feet  and 
^OFFSET  I  (as  described  in  Appendix  D)  is  less  than  15  feet. 

If  these  conditions  are  met,  the  value  of  for  that  cell 

must  be  set  to  one  for  all  values  of  X... 

M 

The  third  iteration  consisted  of  making  further  adjust¬ 
ments  to  the  function  in  an  attempt  to  correct  some  of  its 
deficiencies  as  observed  in  the  plots  in  Appendix  H.  The 
result  of  this  iteration  was  the  following  final  form  of  the 
characterization: 


PK  ’  F  *  e 


X.,  +  a  , 

-(-V— : )6 


(3.5) 


where 


F  -  the  lesser  of  { 


1.22 


.  014R. 


44 


and 

a  and  b  are  determined  as  described  in  Appendices 
D  and  E  with  P  =  .90. 

The  vulnerable  length  is  between  BTP  and  ATP,  where 


ATP  =  5 


and 


CH 


BTP  =  [-9.5  +  2.5(sin[|])  +  2CH] 


is  the  greater  of  t  /12. 5 


and  the  dropout  bias  (see  Appendix  E )  is  1  foot. 


The  13  plots  of  this  final  characterization  correspond¬ 
ing  to  those  in  Appendix  H  are  shown  in  Appendix  I.  Also, 
some  selected  plots  of  this  final  characterization  versus  the 
original  data  are  contained  in  Appendix  J.  The  reader  may 
wish  to  compare  these  to  the  corresponding  plots  in  Appendix  F 


Assessment  of  the  Fit 

The  13  combinations  of  \p ,  VR,  g,  and  o  represented  by 
the  plots  in  Appendix  I  all  have  values  of  ^ ,  VR,  and  B  which 
were  used  to  generate  the  original  data.  The  ultimate  test  of 
this  entire  method  is  to  compare  the  final  characterization  to 
SHAZAM  at  points  in  the  region  of  experimentation  which  were 
not  used  to  generate  the  function.  The  purpose  of  the  test  is 
to  estimate  the  total  error  caused  by  the  use  of  the  character 
zation  in  lieu  of  SHAZAM. 

As  previously  mentioned,  there  is  no  practical  way  to 
precisely  calculate  this  error.  It  is,  therefore,  difficult 
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to  construct  a  test  to  measure  its  value.  The  procedure 
finally  arrived  at  was  to  "randomly"  select  ten  combinations 
of  il» ,  VD ,  and  6  from  the  region  of  experimentation  shown  on 
Table  III  (page  23).  The  value  of  o  was  selected  randomly 
from  the  interval  (5,30).  For  each  of  the  ten  combinations 
of  ip,  VR,  3,  and  o,  six  equi-spaced  values  of  were  exa¬ 
mined  using  SHAZAM.  The  minimum  number  of  iterations  was 
increased  from  30  to  150  and  the  maximum  number  from  70  to 
3S5  in  order  to  reduce  the  halfwidth  of  the  90  percent  confi¬ 
dence  intervals  to  approximately  .04. 

The  ten  selected  combinations  are  shown  in  Table  VI. 

The  reader  should  bear  in  mind  that  a  truly  random  combina¬ 
tion  of  ip,  VR,  and  8  has  a  low  probability  of  representing 
a  trajectory  which  is  achievable  in  reality.  rur  example, 

VR  tends  to  be  higher  for  low  values  of  ip  than  for  high  values 
of  if>.  Therefore,  some  selectivity  was  necessary  in  choosing 
the  combinations.  The  most  important  feature  of  the  selections 
is  that  they  represent  a  cross  -  section  of  the  achievable  tra¬ 
jectories  within  the  region  of  experimentation. 

Trajectory  number  8  was  chosen  for  use  in  an  additional 
experiment.  In  this  experiment,  SHAZAM  was  run  with  angle  of 
attack  (u)  treated  as  a  normally  distributed  random  variable 
with  a  mean  of  5°  and  standard  deviation  of  1°.  The  characteri¬ 
zation  was  then  exercised  also  treating  ot  as  a  random  variable 
with  the  same  distribution.  In  one  case,  a  single  selection 
of  cx  was  made  for  each  cell.  In  the  second  case,  the  "109 
Cell  Model"  was  replaced  by  a  random  sample  of  1000  trajectories 


TABLE  VI 


Ten  Selected  Random  Trajectories 


Traj  No. 

4>° 

VR( fps) 

6° 

o(ft) 

'■‘c 

VM(fps) 

V^,(fps) 

1 

9 

2480 

4 

8.5 

13 

1724 

769 

2 

30 

2604 

14 

7 . 6 

44 

1874 

907 

3 

55 

2447 

24 

16.1 

79 

2042 

1014 

4 

64 

2105 

17 

24.8 

81 

1916 

•  623 

5 

89 

1152 

41 

16.8 

130 

1504 

987 

6 

105 

1799 

21 

13.5 

126 

2148 

797 

7 

122 

1068 

30 

9 . 0 

152 

1929 

1137 

8 

128 

1171 

20 

15.9 

148 

1741 

756 

9 

146 

1467 

11 

29.7 

157 

2099 

716 

i  10 

168 

1646 

5 

15  .  2 

173 

2808 

1177 

Kith  R.^,  6,  and  a  treated  as  random  variables  with  distribu¬ 
tions  as  previously  described. 

The  results  for  the  ten  combinations  are  displayed  in 
the  form  of  plots  in  Appendix  K.  For  the  nine  combinations 
with  c  -  0°  approximately  75  percent,  of  the  characterization’s 
estimates  are  within  .05  of  SHAZAM's  estimates  at  those  points 
where  SHAZAM  was  exercised.  At  only  two  (4£)  of  the  SHAZAM 
estimates  is  the  characterization's  estimate  farther  away 
than  0.1.  The  larger  errors  seem  to  be  a  reflection  of  the 
tendency  of  the  characteri cation  to  drop  from  its  maximum 
earlier  than  the  SHAZAM  estimates. 

For  the  combination  with  a  treated  as  a  random  variable. 
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all  estimates  are  within  -05  of  SHAZAM's  estimates  for  both 
the  "109  Cell  Model"  and  the  random  sample  of  1000  trajec¬ 
tories.  The  only  detectable  difference  between  these  two 
methods  is  in  the  region  between  X^  =  -15  and  =  -5.  The 
maximum  discrepancy  is  about  .02.  Hence,  it  appears  that  no 
significant  error  is  introduced  by  the  use  of  the  "109  Cell 
Model"  with  a  single  random  selection  of  a  for  each  cell. 

Thus,  the  semi-empirical  approach  resulted  in  a  charac¬ 
terization  which  met  the  previously  stated  tolerance  criteria. 
It  remains  to  develop  a  general  procedure  for  characterization 
based  on  this  experience. 
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IV.  Suggested  Method  of  Characterizing 
Endgame  Models  (MOCEM) 


Recall  that  the  objectives  of  this  thesis  are  twofold. 
Eq  (3.5)  and  its  associated  parameters  represents  the  accom¬ 
plishment  of  the  first  objective  of  achieving  a  characteriza¬ 
tion.  The  second  objective,  that  of  defining  a  procedure  for 
arriving  at  a  characterization,  is  discussed  in  this  chapter. 

The  most  desirable  form  of  such  a  procedure  would  be 
that  of  an  "algorithm"  so  that  an  endgame  model  could  be  in¬ 
serted  in  one  side  and  its  characterization  would  then  emerge 
from  the  other  side.  Unfortunately,  the  procedure  used  to 
arrive  at  the  final  characterization  in  Chapter  III  does  not 
easily  lend  itself  to  mechanization.  The  degree  of  subjecti¬ 
vity  used  in  each  iteration  was  so  large  that  it  defies 
rigorous  definition. 

Consequently,  based  on  the  experience  represented  by 
Chapter  111,  only  a  heuristic  procedure  can  be  suggested. 

The  following  steps  are  recommended  to  the  reader  who  wishes 
to  perform  his  or  her  own  characterization  of  an  endgame 
model : 

1)  Disable  the  blast  kill  calculation  portion  of  the 
endgame  program. 

2)  Select  the  set  of  factors  which  seems  to  have  the 
most  direct  impact  on  the  P„  predictions  of  the 
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endgame  program  under  consideration. 

3)  Select  a  range  for  each  parameter  to  define  a 
region  of  experimentation.  Beware  that  the  larger 
the  region  of  experimentation,  the  more  the  prob¬ 
able  effort  required  to  arrive  at  a  characterization. 

4)  Perform  a  small  preliminary  experiment  over  the 
region  of  experimentation  to  determine  the  general 
nature  of  the  function  required.  Select  discrete 
values  of  miss  distance- -i . e. ,  do  not  treat  it  as 
a  random  variable  in  this  step. 

5)  Search  for  a  function  which  exhibits  the  appropri¬ 
ate  behavior.  Hopefully,  a  function  of  the  form: 

V  +  a  7  ji 

f  -  e'1^"'  ,  n  an  integer 

is  satisfactory.  (The  remainder  of  this  procedure 
assumes  that  this  is  the  case.) 

6)  Adjust  the  cone  approximation  of  Appendix  c  to  fit 
the  behavior  of  the  warhead  model. 

7)  Generate  a  data  base  with  the  endgame  program. 

Select  factor  combinations  irom  all  parts  of  the 
region  of  experimentation.  It  is  not  necessary  to 
perform  a  full  factorial  experiment.  Simply  select 
enough  factor  combinations  to  create  what  seems  to 
be  a  representative  sample  of  the  region  of  experi¬ 
mentation.  Again,  do  not  treat  miss  distance  as 

a  random  variab le- - i . e .  ,  select  discrete  relative 
traj  ector ies . 
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8)  Use  the  cone  solutions  of  Appendix  D  (suitably 
generalized  to  account  for  elevation  angle  if  needed) 
to  establish  a  vulnerable  length  for  each  factor 
combination.  Depending  on  the  nature  of  the  warhead 

■  model  and  the  vulnerability  model,  it  may  be  nec¬ 
essary  to  define  more  than  one  vulnerable  length 
along  the  target. 

9)  Observe  the  effects  of  factor  changes  on  the  end¬ 
points  of  the  vulnerable  length(s) .  Use  any  avail¬ 
able  means  to  describe  these  relationships  analyti¬ 
cally.  When  it  is  felt  that  a  sufficient  amount  of 
the  behavior  of  the  vulnerability  model  has  been 
captured  in  this  way,  plot  the  function  versus  the 
data  and  note  any  deficiencies.  Attempt  to  correct 
these  deficiencies  before  proceeding  to  step  10. 

10)  Select  a  "comfortable"  number  of  factor  combinations 
from  those  used  to  generate  the  data.  Using  these 
combinations,  exercise  the  endgame  program,  treating 
TUiss  dist an c c  as  a  random  vari ab le.  be  sure  that 
the  method  of  dealing  with  direct  hit  is  the  same 
for  the  characterization  as  it  is  for  the  endgame 
model.  If  the  distribution  of  miss  distance  is 
circular  normal,  use  the  "109  Cell  Model "  in  Appendix 
G.  Otherwise,  draw  a  sample  of  about  size  300  from 
the  appropriate  distribution  for  each  combination 
of  factprs.  In  either  case,  compare  the  results 
and  attempt  to  correct  any  unacceptable  def iciencies . 
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If  the  deficiencies  cannot  be  corrected,  go  back 
to  step  9;  otherwise  proceed  to  step  11. 

11)  Select  a  number  of  "random'’  combinations  of  factors 
from  the  region  of  experimentation.  Be  sure  to 
include  the  parameter(s)  which  describe  the  miss 
distance  distribution  in  the  random  selection  pro¬ 
cess.  The  number  of  samples  selected  depends  on 
the  subjective  "confidence"  desired.  Exercise  both 
the  endgame  program  and  the  characterization,  treat¬ 
ing  miss  distance  as  a  random  variable.  Again, 
compare  results.  If  unsatisfactoiy ,  go  back  to 
step  9.  Otherwise,  it  is  now  acceptable  to  use  the 
characterization  arrived  at  for  whatever  purpose  it 
was  intended.  (Note:  Any  subsequent  arrival  at 
step  11  will  require  generation  of  a  new  "random" 
sample  of  factor  combinations.  It  is  improper  to 
use  the  same  sample  to  test  for  the  same  deficien¬ 
cies  which  were  discovered  in  prior  accomplishments 
of  step  11.  It  is,  however,  acceptable  to  use  the 
previous  "random"  sample(s)  from  step  11  in  subse¬ 
quent  accomplishments  of  step  10.) 

12)  Be  wary  of  using  the  characterization  outside  of 
the  region  of  experimentation.  Larger  than  accept¬ 
able  errors  may  be  encountered.  It  is  also  prudent 
to  question  any  unusual  behavior  of  the  characteri¬ 
zation  within  the  region  of  experimentation.  Such 
behavior  should  be  verified  with  the  endgame  program 


5  2 


upon  which  the  characterization  is  based. 

There  are  two  special  considerations  to  be  reckoned 
with  in  this  procedure.  The  first  concerns  the  blast  kill 
mechanism.  If  the  experimenter  feels  that  there  are  a  signi¬ 
ficant  number  of  factor  combinations  which  would  result  in 
blast  kills  but  would  not  otherwise  result  in  fragmentation 
kills  in  the  endgame  program,  a  blast  envelope  calculation 
should  be  performed  in  the  characterization  immediately 
following  the  direct  hit  calculation  and  the  blast  kill  calcula¬ 
tion  in  the  endgame  program  should  not  be  disabled.  This  may 
pose  some  difficulty  if  the  actual  physics  of  the  blast  is 
modeled  in  the  endgame  program.  There  should,  however,  be 
no  difficulty  with  the  standard  "envelope"  type  of  calculation. 

The  second  consideration  concerns  the  fuze.  If  a 
fuze  model  is  included  in  the  endgame  program,  it  must  be 
detached  and  dealt  with  separately.  Otherwise,  XM  cannot  be 
treated  as  a  factor  (as  this  method  requires)  since  the  value 
is  determined  by  the  fuze  model.  Upon  completion  of  the 
characterization,  the  fuze  model  can  then  be  linked  to  the 
characterization  in  the  same  manner  that  it  was  originally 
linked  to  the  endgame  program. 


53 


V .  Summary,  Conclu s ions , 
and  Recommendations 

The  objectives  of  this  research  have  been  to:  1) 
characterize  an  endgame  model,  and  2)  develop  a  method  for 
accomplishing  such  a  characterization. 

In  accomplishing  the  first  objective,  a  version  of 
the  endgame  program  known  as  SIIAZAM  was  taken  to  be  repre¬ 
sentative  of  such  programs  in  terms  of  difficulty  in  charac¬ 
terization.  The  factors  studied  were  selected  based  on  the 
geometric  fundamentals  of  the  endgame  problem.  No  experi¬ 
mentation  was  performed  for  the  purpose  of  factor  selection. 

A  purely  empirical  approach  was  pursued  in  the  first 
attempt  to  characterize  SHAZAM .  A  mathematical  function  was 
found  which  demonstrated  the  same  basic  behavior  as  the  data 
generated  with  SHAZAM.  The  high  order  of  interaction  between 
the  factors  prevented  the  accomplishment  of  a  complete  charac 
terization  with  this  approach. 

A  semi- empirica 1  approach  was  then  pursued.  This 
approach  attempted  to  account  for  the  factor  interactions 
caused  by  the  combined  effects  of  the  geometry  and  the  war¬ 
head  characteristics.  It  was  found  that  a  "vulnerable  length 
along  the  target,  centerline  could  be  used  as  a  surrogate  for 
the  vulnerability  model.  Adjustments  to  this  vulnerable 
length  due  to  factor  changes  were  arrived  at  empirically. 
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The  resulting  information  was  used  to  adjust  the  basic  function 
discovered  in  the  empirical  approach.  The  final  characteriza¬ 
tion  resulting  from  the  semi-empirical  approach  was  compared 
to  a  '‘random"  selection  of  data  generated  by  SHAZAM.  The 
characterization's  estimates  were  within  0.1  of  SHAZAM' s  esti¬ 
mates  for  97  percent  of  these  cases. 

The  second  objective  of  developing  a  characterization 
method  was  accomplished  based  on  the  semi-empirical  approach. 
The  resulting  method  was  a  heuristic  one  due  to  the  subjective 
nature  of  the  characterization  process. 

The  amount  of  computer  time  consumed  by  the  semi- 
empirical  approach  in  arriving  at  the  final  characterization 
was  approximately  one  hour  (on  a  CYBER  175)  .  Another  hour 
was  consumed  in  generating  the  60  data  points  used  to  verify 
the  performance  of  the  characterization.  A  computer  program 
implementing  the  mathematics  of  Appendices  C>  D,  E,  and  G 
consisted  of  less  than  150  FORTRAN  statements  and  could 
generate  a  single  Pj,  prediction  in  approximately  0.Z  CP 
seconds.  Since  no  attempt  was  made  to  make  the  program 
efficient,  it  follows  that  even  less  time  is  actually  required. 
Using  the  "109  Cell  Model"  of  Appendix  G,  any  number  of  burst 
points  can  be  evaluated  with  very  little  additional  expendi¬ 
ture  of  computer  time.  For  example,  each  curve  representing 
the  characterization  in  the  plots  in  Appendices  H,  I,  and  K 
consists  of  100  burst  points.  The  total  CP  time  expended  to 
calculate  each  curve  was  0.6  seconds,  including  the  time 
required  to  generate  the  plot.  Thus,  the  time  expended  per 
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I'j.  estimation  has  an  upper  bound  of  0.2  CP  seconds,  but  can  be 
lower  depending  on  the  efficiency  of  the  program  and  the  num¬ 
ber  of  points  calculated. 

Although  it  is  difficult  to  calculate  the  exact  error 
which  results  from  the  use  of  the  characterization  in  lieu  of 
SHAZAM,  the  results  indicate  that  the  error  is  of  approximately 
the  same  magnitude  as  the  uncertainty  in  the  SHAZAM  results 
when  30  to  70  iterations  are  used.  Since  SHAZAM  requires  from 
10  to  20  CP  seconds  for  each  such  P^  estimate,  the  time  savings 
resulting  from  the  substitution  of  the  characterization  for 
SHAZAM  is  approximately  two  orders  of  magnitude  and  possibly 
more . 

Not  a  great  deal  of  time  was  spent  in  attempting  to 
capture  the  details  of  the  vulnerability  model's  behavior. 

It  is  conceivable  that  if  one  was  willing  to  spend  the  time, 
sufficient  refinements  could  be  made  to  the  characterization 
to  virtually  duplicate  the  SHAZAM  results  with  no  detectable 
error . 

Hence,  what  has  been  shown  is  that  an  endgame  model 
can  be  sufficiently  characterized  under  somewhat  limited  con¬ 
ditions.  Whether  every  endgame  model  can  be  characterized  in 
this  way  is  another  question.  Differences  in  warhead  charac¬ 
teristics,  target  characteristics,  arid  endgame  program  metho¬ 
dology  may  have  a  critical  impact  upon  the  ability  of  the 
method  to  accomplish  an  adequate  characterization.  However, 
the  suggested  procedure,  having  worked  in  this  case,  should 
be  the  first  approach  taken  in  any  future  efforts  of  this 
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nature . 


There  are  several  potential  areas  of  future  research 
which  would  provide  more  generality  to  the  results  presented 
herein.  The  inclusion  of  elevation  angle  is  certainly  the 
area  which  demands  first  attention.  In  addition,  application 
of  the  method  to  different  combinations  of  endgame  models, 
warhead  models,  and  vulnerability  models  would  illuminate  the 
question  of  general  applicability  of  the  approach.  Also,  a 
method  for  mechanizing  the  suggested  characterization  proce¬ 
dure  is  highly  desirable.  Accomplishment  of  such  an  algorithm 
depends  highly  on  the  development  of  a  "measure-of-merit”  of 
fit  which  can  suitably  replace  the  subjective  judgment  as  used 
herein.  Finally,  expansion  of  the  method  of  characterization 
to  include  the  fuze  model  should  be  investigated. 
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APPENDIX  A 

Relationship  Between  Selected  Factors 
and  S HA 2AM  Inputs 


This  appendix  contains  a  description  of  the  relation¬ 
ship  between  the  selected  factors  VR,  6,  and  ij;  to  the  SHAZAM 
inputs  V„,  Vy,  ^  when  X  =  0. 

From  Figure  A-l,  and  if  a  =  0,  we  have: 


=  ip  +  6 


Given  VD  ,  B ,  and  4),  by  the  Law  of  Sines, 
K 


Similarly , 


V 


M 


VT 


sin  ij;  _  sin(180-^c) 


=>  V„  = 


VR  sin  ^ 


M  sm 


Vp  sin  B 
v  =  — . - 

T  sm  Vr 


Given  VT,  VM,  and  ipc ,  by  the  Law  of  Cosines, 


V 


V  *  VT2  ■  ZVMVT  =05(180  -  *c) 


R 


V  =  sin  '1  [(VM/VR)  sin  i|»c] 


8  =  4’c  ~  i|> 


If  a  t  0,  B -cx  should  be  substituted  for  B  in  the  above  expres¬ 
sions.  The  values  of  X^,  1^,  and  Z^  are  the  same  for  the 
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selected  factors  as  for  SHAZAM.  and  6  are  defined 
Chapter  III. 


APPENDIX  B 


Results  of  the  Empirical  Approach 

This  appendix  contains  plots  which  depict  the  results 
of  the  empirical  approach  discussed  in  Chapter  III, 

Figures  B-l  through  B--5  represent  the  plots  of  "a" 
versus  e  (location  function)  for  the  five  values  of  in  the 
first  stage  of  the  experiment.  Figures  B~6  through  B-10  are 
the  plots  of  "b"  versus  6  (halfwidth  function)  for  the  five 
values  of  iJj.  The  curves  are  the  "fits"  whose  equations  are 
shown  on  the  individual  plots. 

The  fit  equations  in  Figures  B-l  through  B-10  are  all 
of  the  form: 

Cx  +  C2  sin2(ie  +  C3) 

The  data  points  in  Figures  B-ll  through  B-16  represent  the 
values  of  C^,  and  as  a  function  of  ^  for  both  the  loca¬ 
tion  function  and  the  halfwidth  function.  is  the  constant 

term,  is  the  amplitude  term,  and  is  the  phase  term. 

The  curves  represent  the  fitted  equation  which  is  shown  on 
each  plot. 

Figures  B-17  through  B-32  are  plots  of  the  characteriza¬ 
tion  (from  the  empirical  approach)  versus  the  data  generated 
by  SHAZAM  for  stage  1  of  the  experiment.  The  values  of  the 
factors  are  shown  on  each  plot.  The  "error  bars"  on  the  data 
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points  represent  the  90  percent  confidence  interval  about  the 
mean  P^- - i. e . ,  ± . 10 . 
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APPENDIX  C 


Approximation  of  the  Warhead  Pattern 
With  a  Cone 

This  appendix  describes  the  cone  approximation  of  the 
warhead  pattern  discussed  in  Chapter  III. 

Let  Vr  be  the  fragment  velocity  and  <J>  be  the  angle  of 

r 

fragment  ejection  as  measured  from  the  missile  centerline 

under  static  conditions.  Assume  that  X  =  o  and  a  =  0.  Then 

in  the  relative  reference  frame,  the  situation  shown  in 

Figure  C-l  represents  the  vector  geometry  in  a  plane  which  is 

parallel  to  the  X-Y  plane  and  contains  the  missile  centerline. 

and  r 2  are  the  angles  (in  this  plane)  between  VR  ana  VR 

and  VD  and  \\  ,  respectively,  and  are  given  by: 

K  h  2 

=  4>  +  B 
n2  =  <J>  -  6 

The  magnitudes  of  Vc  and  are  both  equal  to  V,. .  V_  repre- 

hl  2  1 

sents  the  velocity  vector  of  the  fragment  which  is  ejected  to 

the  port  side  of  the  missile  parallel  to  the  X-Y  plane.  VP 

2 

represents  the  analogous  vector  on  the  starboard  side. 

As  shown  in  Figure  C-2,  the  vector  addition  of  and 
VP  (either  VP  or  VP  )  yields  a  resultant  VP’  which  makes  an 
angle  V  with  VR  where 
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VpSin  n 

u  =  tan  (VR+VFco^ 

Figure  C-3  depicts  the  relationship  between  these 

resultants,  Vn ,  and  the  missile  centerline.  The  cone  half- 
K 

angle,  6,  is  given  by: 


The  angle  e  between  the  cone  axis  and  VR  is  given  by: 

e  =  6  -  y2 

Therefore,  the  angle  between  the  cone  axis  and  the  mis¬ 
sile  centerline  is  just  6  -  e.  If  a  is  not  zero,  then  3  +  a 
is  substituted  for  S  in  all  expressions. 

Figures  C-4  through  C-8  are  comparisons  between  this 
cone  approximation  and  the  true  dynamic  warhead  pattern.  The 
plane  of  the  plot  represents  the  Y  -  Z  plane  in  the  target 
coordinate  system.  The  axes  are  centered  on  the  target  centroid. 
The  squares  represent  the  penetration  points  of  the  plane  for 
the  true  dynamic  warhead  pattern  if  fragments  were  spaced 
every  10°  around  the  warhead.  The  +'s  are  the  corresponding 
points  of  intersection  between  the  approximating  cone  and  the 
Y-Z  plane  of  the  target  coordinate  system.  The  scales  are  in 
units  of  feet,  and  DP  is  the  distance  in  feet,  measured  along 
the  missile  centerline,  to  the  Y-Z  plane.  For  each  of  these 
plots,  Vp  =  5500  feet  per  second  and  $  =  84°.  The  values  of 
the  other  factors  are  shown  on  each  plot. 
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The  reader  should  remember  that  it  is  the  continuous 
intersection  with  the  Y-Z  plane  which  is  of  interest,  not  the 
individual  points  of  intersection.  Therefore,  that  the  appro¬ 
ximation  calculates  individual  fragment  penetrations  of  the 
plane  incorrectly  is  not  important- -  it  is  the  pattern  which 
is  of  concern  in  the  semi-empirical  approach  of  Chapter  III. 


Fig  C-l:  Vector  Geometry 
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APPENDIX  D 


Simultaneous  Solution  of  the 
Cone/Line  Intersection 

This  appendix  develops  the  mathematics  required  to 
solve  for  the  intersection  between  a  line  and  a  cone.  In 
the  context  of  this  thesis,  the  line  represents  the  infinite 
line  which  is  coincident  with  the  target  centerline  and  the 
cone  represents  the  warhead  pattern  approximation  presented 
in  Appendix  C.  The  value  of  X  is  assumed  to  be  zero  for  all 
of  the  derivations  contained  herein. 

If  ZM  =  0,  Figure  D-l  represents  the  geometry  of  the 
problem.  There  are  two  ways  of  solving  for  the  intersection, 
depending  upon  the  information  given.  If  DT  is  known,  then 
it  is  possible  to  solve  for  the  distances  D1  and  D2.  Con¬ 
versely,  if  either  D1  or  D2  is  known,  it  is  possible  to  solve 
for  DT .  In  either  case,  the  following  relationships  hold: 

VM>  hi  are  as  dcfincd  in  Chapter  III 

XM  -  -(DT  t  DF) 

DF  =  Y^/tan'J' 
w  =  -90  +  (p  +  c) 
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^OFFSET  ~  Uni’ 

Define  the  primed  coordinate  system  to  have  its  origin 
at  the  warhead  centroid,  its  X-axis  (X')  coincident  with  the 
cone  axis,  its  Z-axis  (Z')  in  the  same  direction  as  the  Z-axis 
of  the  target  coordinate  system,  and  its  Y-axis  (Y*)  defined 
according  to  the  right  hand  rule.  Then,  if  the  coordinates 
of  a  point  on  the  X-axis  in  the  target  coordinate  system  are 
(X,  0,  0),  then  its  coordinates  in  the  primed  system  are  given 

by: 

X'  =  (X  -  A)  sin  u  +  DP 
Y '  =  (X-A)  cos  w 

7  i  —  .  7 

"M 

where 

A  =  n3  +  DOFFSET 

In  the  primed  coordinate  system,  the  equation  of  the  warhead 
cone  is 

Y’Z  +  Z,Z  =  X’z  (tan  6)" 

To  solve  for  D1  and  D2  when  DT  is  known,  let  (X1 ,  0,  0) 
and  C  X  ^ »  0,  0)  be  two  arbitrary  points  on  the  X-axis  of  the 
target  coordinate  system.  Let  (X^',  Y^',  Z^')  and  (X^',  Y 2 ’  » 

Z , 1 )  be  the  coordinates  (respectively)  of  these  two  points  in 
the  primed  coordinate  system.  Then,  if  TP  is  any  point  on  the 
X-axis  of  the  target  coordinate  system,  there  exists  a  t  such 
that  the  coordinates  of  TP  in  the  primed  coordinate  system  are 
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If  TP  lies  on  the  cone,  then 

ytf'2  +  2tp'2  =  xXP,2Ctan  6)2 

=>  (Y  x  '  -  ( V  2  ’  -  Y  x  ’  )  t )  2  +  (-ZM)2  =  CX1’  +  a2’-X1')t)2tan26 
=>  Y1’2  +  2Y1*(Y2--Y1’)t  -  (Y2,-Y1,)2t2  +  Z^2 

=  tan26  (X1'2  +  2X,’ (X2’-X1,)t  +  (X7  '  -  X  '  ft2) 

After  some  rearrangement,  this  yields: 

[(Y2’-Y1*)2  -  (X2'-X1')2tan26)t2 

+  2  [Y1 ’ CY2'-Y1')  -  X1'(X2’-X1')tan26]t 

+  [ZM2  +  Y ^ ' 2  -  X1'2tan26]  =  0 
which  has  solutions: 


■B±  / B  -  4AC 


where 


-  f  v  i  „  v  i  2 

-  1 1  2  >  i  1 


7  7  , 

v  -V  J  )  l  Clll  w 


B  «  2  [  Y  x ' ( Y  2 '  -  Y1')  -  X1'tan26(X2’  -X1 1 ) ] 

7  7  2  7. 

C  =  Z„  +  Y, ’  -  X, ’  tanz6 

Ml  1 

2 

If  B  <  4AC ,  the  cone  and  the  X-axis  of  the  target  coordinate 

2 

system  do  not  intersect.  If  B  =  4AC,  then  the  X-axis  of  the 

target  coordinate  system' is  tangent  to  the  cone.  Otherwise, 
n 

B “  >  4AC  and  there  are  two  solutions:  t,  and  t„. 


The  value  of  t  is  just  the  parameter  associated  with 
the  parametric  description  of  the  X-axis  of  the  target  coordin¬ 
ate  system.  If  the  two  points  defining  the  line  are  the  same 
in  absolute  space,  the  value  of  t  is  the  same  in  any  coordinate 
system.  Since  (X^',  Y^',  Z  n  ' )  and  (X9',  Y/,  Z^')  correspond 
to  X^  and  X9  ,  respectively,  in  the  target  coordinate  system, 
then  the  coordinates  of  ATP  and  BTP  in  the  target  coordinate 
system  are  just 

(X^TP ’  0 >  and  ^ 

where 

XATP  X1  (-X2'Xl/)tl  and  XBTP  =  X1  +  (X2'X1^2 

There  is  a  pathology  associated  with  this  solution.  The 
mathematical  equation  of  a  cone  actually  describes  two  cones 
in  mathematical  space.  The  warhead  cone  is  that  portion  of 
the  mathematical  cone  for  which  the  X’  coordinate  is  positive. 
Therefore,  any  of  the  solutions  for  t  which  result  in  a  nega¬ 
tive  •'ralue  of  X’  should  be  rejected.  Note  that  in  the  situa¬ 
tion  depicted  in  Figure  D-l,  the  intersection  of  the  cone  and 
the  X-Z  plane  of  the  target  coordinate  system  is  an  ellipse. 

If,  and  only  if,  this  is  the  case  will  both  ATP  and  BTP  exist 
for  a  given  value  of  DT. 

Note  that  in  Figure  D-l,  =  0.  If  Z^  ^  0,  then  the 
lines  connecting  ATr  and  BTP  to  the  origin  of  the  primed  system 
will  not  make  an  angle  5  with  the  cone  axis  in  the  X'  -  Y'  plane. 
Rather,  the  total  angle  (in  three  dimensions)  between  the  cone 
axis  and  either  of  the  two  lines  is  6,  since  by  assumption. 
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ATP  and  BTP  lie  on  the  cone. 

Now,  if  ATP  and  BTP  are  specified,  it  is  possible  to 
solve  for  DT.  For  purposes  of  generality,  assume  that  ATP 
and  BTP  need  not  lie  on  the  cone  s imultaneously- - i . e .  for  a 
given  value  of  DT.  The  following  development  centers  on  the 
solution  for  DT  if  ATP  is  known.  The  development  if  BTP  is 
known  is  completely  analogous. 

Let  the  coordinates  of  ATP  in  the  target  coordinate 
system  be  (xAyp >  0»  0) ,  then  its  coordinates  in  the  primed 
coordinate  system  are: 

X*  =  ( XATP  '  ^  sin  to  +  DP 


( Y  _  DT  sin  e  _  M  ■,  ■  .  DT  sin  \p 

ATP  sin (iii  +  e)  s in(w)  '  w  T  sin(\I;  +  e) 

„  v  sin  “  j  -  sinesinw^ 

-  XATpsinU  -  -sa-f-  *  DT[  ^ - ) 


It  can  be  shown  that 


sin  jp  -  sin  e  sin  ui 
sin(^  +  e) 


cos  E 


Thus,  the  equation  for  X'  becomes 
YMsinw 

■  XATP  Sin  "  •  “sirnj;  *  DT  cos  E 


y '  -  cx, 


A)  cos  w 


=  (XATP  '  sin  O  +  e) 

rv  DT  sine  YM  . 

(XATP  '  sinTlFH'  _  sin^  5in^  T 


=  X^p  sin(ijj+e)  -  DT  sin  c  - 


Ym  sin  (iJj  +  e) 
sin  \{j 
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Now,  substituting  these  expressions  into  the  cone  equation, 


where 

2  2  2 
A  =  sin  £  -  tan  <5  cos  e 

7 

B  =  -  2 (X.Tp-YM/sin^) [sin£  cosw  +  tan  6  sinw  cose  ] 
C  =  Z^2  +  (X^pp- Y^/sinifO  2  [  cos  2w  -  tan2<5  sin2w] 


Again,  there  are  pathologies  associated  with  this  solution 

when  B2  ^  4AC .  If  B2  <  4AC ,  there  is  no  value  of  DT  for  which 

2 

ATP  lies  on  the  cone.  If  B  =  4AC ,  there  is  only  one  value  of 
DT  for  which  ATP  lies  on  the  cone.  Otherwise,  there  are  two 
values  of  DT  for  which  ATP  lies  on  the  cone.  Since 


XM  =  -  DT  +  DF 

choose  that  value  of  DT  which  minimizes  X...  This  assures  that 

M 

the  forward  half  of  the  mathematical  cone  (i.e.  the  warhead 
cone)  is  the  one  on  which  ATP  lies. 
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Fiy  D-l:  Geometry  of  the  Cone/Line  Intersection  Problem 


APPENDIX  E 

Determination  of  "a"  and  "b1 


This  appendix  describes  how  a  and  b  are  determined 
given  the  values  of  X^  at  which  ATP  and  BTP  are  intersected 
by  the  warhead  cone. 

Let  X..  be  the  value  of  X,.  when  ATP  lies  on  the  leading 

mf1  m 

warhead  cone,  be  the  value  of  X^  when  ATP  lies  on  the  trail¬ 
ing  warhead  cone,  X^  be  the  value  of  when  BTP  lies  on  the 

leading  warhead  cone,  and  X.„  be  the  value  of  XM  when  BTP  lies 

KR2  M 

on  the  trailing  warhead  cone.  Xw  ,  X..  ,  XM  ,  and  XM  are 

MF1  mF2  R1  R2 

determined  as  described  in  Appendix  D.  In  this  context,  ATP 

refers  to  the  forward  limit  of  the  vulnerable  length  (i.e.  the 

end  of  the  vulnerable  length  closest  to  the  nose  of  the  target) 

while  BTP  refers  to  the  rearward  limit  (see  Figure  E-l).  Let 

P  denote  the  proportion  of  maximum  P^  which  is  achieved  when 

either  ATP  or  BTP  is  first  or  last  contacted  by  a  warhead  cone. 

Define  XM  and  XM  as  follows: 

‘  F  MR 


if  XM  >  XM 

mfi  mri 


Otherwise , 


XM  =  XM 
Mp  Mp2 


an<1  \  ■  \ 


Up  Mpi 
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This  definition  assumes  that  the  value  of  first 
reaches  P  •  P,.  whenever  the  vulnerable  length  is  first  con- 

kmax 

tacted  by  a  warhead  cone,  whether  it  be  on  the  ATP  end  or  the 
BTP  end.  P„  is  then  last  equal  to  P  •  Pv  whenever  the  vul- 

K  kmax 

nerable  length  is  last  contacted  by  a  warhead  cone,  again 
regardless  of  which  end  (see  Fig.  E-2). 


Then , 


Now , 


Let  W  =  |X^  -  X^  J  and  L  =  the  minimum  of  X^  and 


a  =  -(L  +  W/2) 


X.+a  A 
.  r  M  ,6 


=  P  whenever  X..  =  X.,  or  X..  -  X.. 

M  Mp  M 


In  particular, 


X,,  Ta 

Mr;  A 
r  b  1 6 

.  r - J 


Va 

->  -  InP  -  (—I—)6 


I 


X, .  ^  a 

f.  _  M-r 

=  >  V-lnP  =  — l — 


=  >  b 


X.,  +  a 

3 _ 

s/TnP 


There  is  a  phenomenon  whereby  the  vulnerable  length  is  some¬ 
times  last  contacted  by  a  warhead  cone  after  both  ATP  and  BTP 
have  last  been  contacted.  This  situation,  depicted  in  Figure 
E-3,  can  occur  only  when  the  intersections  of  the  warhead 
cones  with  the  X-Z  plane  in  the  target  coordinate  system  are 
ellipses.  The  situation  is  of  concern  whenever  a  point  bet¬ 
ween  ATP  and  BTP  becomes  a  point  of  tangency  between  the  line 


.  K 


r  v 


and  a  warhead  cone.  The  result  is  that  the  vulnerable  length 
suddenly  "drops- out"  of  the  fragment  pattern. 

Mathematically,  the  necessary  condition  for  dropout  is 
that  B  =  4AC  where  A,  B,  and  C  are  as  defined  by  the  first 
solution  in  Appendix  D.  Substitution  of  these  expressions 
for  A,  B,  and  C  into  B  =  4AC  and  subsequent  algebraic  mani¬ 
pulation  reveals  that  an  equivalent  necessary  condition  is 
that 

DT  =  |  Z^/sinip  |  /(cos ^Wtan^6) -sin^io 


Thus,  to  account  for  drop-out,  one  should  substitute 


this  value  of  DT  (say  DT 


DROPOUT 


)  into  the  first  cone/line 


intersection  solution  of  Appendix  D.  If  and  only  if  the 

resulting  point  is  between  BTP  and  ATP  does  drop-out  occur. 

r f  xt  does’  then  XM  =  dtdropout  and  Xw  the  minimum  of 

K  F 

9  X» .  %  and  XCw  • 

MF1  HFZ  MR1  HR2 

When  the  vulnerable  length  drops-out  of  the  warhead 
cone,  the  fragment  pattern  is  still  intersecting  the  target 
fuselage  surrounding  it.  Obviously,  there  are  vulnerable  com¬ 
ponents  in  this  part  of  the  fuselage,  for  otherwise  the  vul¬ 
nerable  length  would  be  elsewhere  along  the  target  centerline. 
Therefore,  it  is  necessary  to  "bias"  X^  by  some  amount  to 
account  for  the  width  of  the  fuselage.  This  is  accomplished 
by  adding  a  "drop-out  bias"  to  E)PDRqP0UT  in  obtaining  X^ 

R 

(See  Fig .  E-4) , 
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Target 

Centroid 


ATP 


^ - 0 


To  Nose 


L:  The  Vulnerable  Length  1 

j 


1 


Leading  Trailing 
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APPENDIX  F 

The  First  Characterization  Applied 
to  the  Basic  Data 

The  48  plots  in  this  appendix  represent  a  comparison 
between  the  first  semi- empirical  characterization  and  the 
basic  data.  The  following  is  a  description  of  this  first 
characterization 


PK  -  F  .  a 


where 


the  lesser  of  { 


1.0 
1 .22 


.  14R 


M 


and  a  and  b  are  determined  as  described  in 
Appendices  D  and  E  with  P  =  0.9  and 
ATP  =  5 

BTP  =  -7+  2.5  CH 


where  CH  is  the  greater  of  { 


0 

(RM- 15) /I 2 . 5 


and  the  drop-out  bias  is  one  foot. 

There  is  one  plot  for  each  entry  in  Table  V.  The  combination 
of  factor  values  is  given  on  each  plot.  The  error  bars  are 
±  .10  90%  confidence  intervals  about  the  SI-1AZAM  estimates. 
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APPENDIX  G 


"109  Cell  Model"  of  the  Circular 
Normal  Distribution 

The  "109  Cell  Model  of  the  Circular  Normal  Distribution 
consists  of  a  series  of  rings  which  are  sub-divided  into  cells 
as  depicted  in  Figure  G-l.  The  probability  of  occurrence  of 
each  cell  in  the  first  seven  rings  is  .01.  The  cells  in  the 
outer  ring  each  have  a  .0025  probability  of  occurrence.  The 
point  in  the  center  of  each  cell  is  the  probabilistic  "center" 
of  that  cell--i.e.  the  probability,  within  that  cell,  of  being 
either  closer  or  farther  from  the  origin  is  the  same  and  there 
is  also  an  equal  probability  of  being  clockwise  or  counter¬ 
clockwise  of  this  point.  The  radius  of  each  ring  in  terms 
of  o  is  given  below. 


Ring 

if  of  Cells 
in  Ring 

Total 

Inside  Ring 

Ring 

Radius ( o) 

Probab  ilistic 
Center  Radius 
(°) 

1 

1 

1 

.  1418 

0 

2 

7 

8 

.  4083 

.  2949 

3 

13 

21 

.6867 

,5556  i 

d 

18 

39 

.9943 

.  8433 

5 

21 

60 

1.3540 

1.1715 

6 

21 

81 

1.8226 

1.5700 

1 

16 

97 

2.643 

2 . 1390 

8 

17. 

109 

CO 

2.9871 
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APPENDIX  II 


The  First  Ch  '.racter  ization  With 
Random  Miss  Distance 

The  thirteen  plots  in  this  appendix  are  a  comparison 
between  data  generated  by  SHAZAM  (the  points)  and  the  cor¬ 
responding  data  (the  curves)  produced  by  the  first  semi- 
empirical  characterization  described  in  Appendix  F.  For  all 
of  these  plots,  miss  distance  was  tieated  as  a  random  vari¬ 
able  in  both  SHAZAM  and  the  characterization.  Values  of  y, 

,  3,  and  o  are  given  on  each  plot. 
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APPENDIX  J 


Selected  Plots  of  the  Improved 
Characterization  Applied  to  the 
Basic  Data 

The  10  plots  in  this  appendix  are  comparisons  of  the 
final  characterization  with  selected  factor  combinations  from 
the  basic  data.  The  corresponding  plots  for  the  first  charac¬ 
terization  appear  in  Appendix  F. 
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to  the  10  Selected  Random  Trajectories 


The  11  plots  in  this  appendix  show  the  final  characteri 
zation  versus  the  data  generated  by  SHAZAM  for  the  ten  random 
trajectories  in  Table  VI.  For  all  of  these  plots,  the  error 
bars  on  the  data  points  are  intended  to  graphically  depict  a 
distance  of  +0.1  from  each  point.  They  do  not  represent 
confidence  intervals  in  this  appendix.  The  90%  confidence 
intervals  about  the  SHAZAM  P^'s  for  these  ten  random  trajec¬ 
tories  is  +0.04. 

Figures  K-8  and  K-9  are  the  trajectory  with  a  treated 
as  a  random  variable.  In  both  plots,  the  SHAZAM  data  is  the 
same.  The  characterization  curve  in  Figure  K-8  was  generated 
with  the  109  Cell  Model  while  in  Figure  9,  it  was  generated 
with  1000  randomly  selected  combinations  of  RM ,  O,  and  a. 
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required  in  any  given  endgame  analysis.  The  main  objective  of 
this  research  was  to  develop  a  method  for  characterizing  an  end¬ 
game  model  using  only  a  fraction  of  the  amount  of  data  generated  _ 

DD  ,  ™Tn  1473  EDITION  OF  1  NOV  65  IS  OBSOLETE  UNCLASS  IFIED 


4f 


SLCUPITY  CLASSIFICATION  OF  THIS  PAGE  fWton  Dei#  E nfaiad) 


_ UNCLASSIFIED _ 

security  CLASSIFICATION  of  this  PAGEfWiAn  Q*ta  Snund) 


^|3Iock  20  : 

by  such  analyses.  In  pursuit  of  this  objective,  a  second  objec¬ 
tive  of  actually  accomplishing  a  characterization  o  an  existing 
endgame  program  was  established. 

."The  endgame  program  SHAZAMwas  used  to  generate  data  upon 
which  to  base  a  characterization.  A  purely  empirical  approach 
was  first  attempted  wherein  all  main  and  interactive  effects  of 
the  critical  factors  were  to  be  assessed  strictly  by  experimen¬ 
tation  with  and  observation  of  the  behavior  of  SHAZAM.  The 
high  order  of  interactions  between  the  factors  precluded  com¬ 
pletion  of  this  approach. 

A  semi -empirical  approach  was  then  pursued  whereby  the 
interactions  between  the  endgame  geometry  and  warhead  charac¬ 
teristics  were  accounted  for  geometrically.  A  two-parameter 
exponsntia.1  function  was  used  to  describe  the  behavior  of  P,, 
as  a  function  of  the  warhead  burst  point.  Values  of  the 
two  parameters  were  calculated  based  upon  the  settings  of  the 
remaining  factors.  A  "vulnerable  length"  was  used  to  account 
for  the  target's  vulnerability.  After  several  iterations  of 
refinement,  the  characterization  was  compared  to  SHAZAM  for  ten 
randomly  selected  trajectories.  The  characterization's  P„. 
estimates  were  within  0.1  of  SHAZAM' s  estimates  for  97 
percent  of  these  random  data.  A  simple  computer  program  imple¬ 
menting  the  mathematics  of  the  characterization  required  a 
maximum  of  0.2  seconds  of  computer  time  per  estimate. 

Approximately  two  hours  of  computer  time  were  consumed  by  SHAZAM 
in  generating  the  necessary  data.  Thus,  with  a  relatively 
modest  investment  of  computer  time,  two  orders  of  magnitude 
improvement  in  time  required  per  ?K  estimate  was  achieved. 

— — ■ — 7  The  suggested  method  for  characterizing  air-to-air  missile 
endgame  programs  is  based  on  the  semi-empirical  approach  and  its 
application  to  SHAZAM.  The  method  is  heuristic,  rather  than 
algorithmic,  due  to  the  subjective  nature  of  the  characteriza¬ 
tion  process . 
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